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ABSTRACT 


A noise model is formulated to describe the impulse noise in many 
digital systems. A simplified model, which assumes that each noise 
burst contains a randomly weighted version of the same basic waveform, 
is used to derive the performance equations for a correlation receiver. 
The expected number of bit errors per noise burst is expressed as a 
function of the average signal energy, signal-set correlation coef- 
ficient, bit time, noise-weighting-factor variance and probability 
density function, and a time range function which depends on the cross- 
correlation of the signal-set basis functions and the noise waveform. 

A procedure is established for extending the results for the simplified 
noise model to the general model. Unlike the performance results for 
Gaussian noise, it is shown that for impulse noise the error perforr^ 
mance is affected by the choice of signal-set basis functions and that 
Orthogonal signaling is not equivalent to On-Off signaling with the 
same average energy. The error performance equations are upper bounded 
by a procedure which removes the dependence on the noise waveform and 
the signal-set basis functions. A looser upper bound which is also in- 
dependent of the noise-weighting-factor probability density function is 
obtained by applying Chebyshev’s inequality to the previous bound result. 
It is shox\m that the correlation receiver is not optimum for impulse 
noise because the signal and noise are correlated and the error perfor- 
mance can be improved by' inserting a properly specified nonlinear device 
prior to the receiver input. 
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CHAPTER 1 


INTRODUCTION 


!•! STATEMENT OF PROBLEM 

Traditionally, the performance analyses of digital systems have 
incorporated additive white Gaussian (AWG) noise to model the inter- 
ference introduced by the communication channel. With this assump- 
tion it has been shown that the optimum receiver consists of a bank 
of correlators, each of which compares the incoming noisy waveform 
with one of the possible received noise-free signals. Although the 
AWG noise assumption gives satisfactory error rate predictions for 
many situations, it has failed to predict the much higher error 
rates observed on some classes of digital systems. Another type of 
additive noise, commonly referred to as impulse noise, is currently 
being modeled and used to explain this discrepancy on atmospheric 
radio channels and both switched and dedicated baseband channels. 

The impulse noise models proposed in the literature differ in sev- 
eral respects, but all impulse noise is characterized by brief 
periods of large amplitude excursions, separated by intervals of 
quiescent conditions . 

The relatively new area of impulse noise analysis suffers from 
the problems of nonstandardization and incompleteness. It is diffi- 
cult to compare the results of two authors and impossible to obtain 
comparison of the major binary signaling techniques, because of the 
multitude of noise models, signal sets, analytical methods and 


1 



2 


simplifying assumptions. In this work the primary objective is to 
establish a mathematical model that can be used to approximate addi- 
tive impulse noise and then apply this model in the error perfor- 
mance analysis for a correlation receiver. The secondary goal is to 
present a collection of specific performance results from which it 
will be possible to draw some general conclusions regarding impulse 
noise and its effect on digital systems. 

1.2 IMPULSE NOISE LITERATURE REVIEW 

Almost all papers concerning impulse noise have been published 
since 1960. A wide variety of mathematical models have been proposed 
to describe different aspects of impulse noise. Only a few authors 
[1,2]* have discussed the noise waveforms obtained from actual systems, 
while much emphasis has been placed on empirically obtained amplitude 
and/or time distributions for the noise [l,3-^8]. Proposed impulse noise 
models have assumed such waveforms as the ideal impulse function, with 
either constant weight [9] or probabilistic weighting [10,11], and the 
unit-impulse response of the channel being studied [12,13]. Another 
approach has been to assume various probabilistic descriptions for the 
noise at the receiver output [14]. The majority of performance studies 
have been limited to the analysis of a specific binary signaling tech- 
nique [2,6,12,13,15-20] although several authors [9,11,21-24] have ob- 
tained performance curves and compared more than one signal set based 


^Enclosed numbers refer to corresponding entries found in List of 
References . 
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on their assumed noise models, Zeimer [10] and Millard & Kurz [25] 
have extended their analyses to the consideration of M-ary signaling. 

The basic approach to error analysis for the practical case of a 
general noise waveform which can change from one noise burst to another 
has been developed by Houts & Moore [11] and forms the foundation of 
the work documented here. Although specific performance results were 
not obtained for the general model, upper and lower bounds on signal 
set error performances were derived. 

Signals corrupted by impulse noise are often detected by receiv- 
ers that were designed for use in AWG noise. It is reasonable to ex- 
pect that the error performance of these receivers can be improved. 
Several receiver-improvement modifications have been proposed [26-30] 
based upon intuition or experim.ental results. The subject of optimum 
receiver design for binary signals in the presence of impulse noise 
has received limited treatment in the literature. Snyder [31] has 
investigated the design of an optimum receiver for VLF atmospheric 
noise. The noise model has a combination of AWG noise and impulse 
noise generated by exciting a filter with idealized delta functions. 
Specific results are not obtained and remain a topic for further 
investigation. A different approach has been used by Rappaport & 

Kurz [32] who assume that the noise waveform duration is short com- 
pared to the bit time and that independent samples can be obtained 


by uniformly sampling the received signal. 
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1 . 3 APPROACH OUTLINE 

The first step in any error performance analysis should be the 
development of a realistic noise model. Simplifying assumptions are 
sometimes included in the model for the purpose of making the analysis 
mathematically tractable, which in turn aids in understanding the role 
played by the various noise parameters. However, the results from such a 
model often fail to predict actual system performance. In Chapter 2 a 
generalized impulse noise (GIN) model is presented which can be used to 
approximate the impulse nois^ structure in a broad class of digital sys- 
tems. In particular, the model is applicable to a baseband channel with 
near dc response of the type found in telephone circuits and sophisti-- 
cated avionics systems such as the proposed data bus for NASA’s Space 
Shuttle vehicle. The performance analysis for a band-limited channel 
using a correlation receiver in the presence of impulse noise described 
by the GIN model is presented in Chapter 3. In addition, new and more 
definitive performance bounds are obtained for typical binary signal 
sets. The effects of selecting various noise weighting factor descrip- 
tions, signal sets, and noise waveforms are documented in Chapter 4. 
Several performance improvement methods are surveyed in Chapter 5 and 
specific results are presented for the case where a nonlinear device 
is inserted prior to the correlation receiver. The results show that 
the performance can be improved; consequently, the correlation re- 
ceiver is not optimum for impulse noise corrupted signals. 



CHAPTER 2 


IMPULSE NOISE MODEL 


A mathematical model is developed which is sufficiently flexible 
to approximate the additive impulse noise in a broad class of digital 
systems. Also, the power-density spectrum and probability density 
function (p.d.f.) are established for certain special cases of the 
noise model. 

2 . 1 HEURISTIC CONSIDERATIONS 

The waveform description of impulse noise depends on the noise 
source, the coupling media, and the characteristics of the channel 
under observation. Capacitive and inductive coupling along with 
switching transients comprise the major methods by which impulse 
noise is introduced into a system. For instance, impulse noise in 
a channel might be caused by a lightning discharge or engine igni- 
tion noise. Telephone and telemetry systems experience noise from 
electromechanical relay switching. Other sources include environ- 
mental conditions such as shock, vibration, and maintenance disrup- 
tions, i.e., a temporary short circuit. These sources have the com- 
mon characteristic that the noise produced occurs in bursts of large 
amplitude variations separated by relatively long intervals of neg- 
ligible amplitude variations. This observation is the foundation 
for the mathematical model presented in Section 2.2. 
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The determination of the noise waveforms can be a time consuming 
and expensive process. The number of different waveforms encountered 
can be very large, and normally the measuring equipment must have 
large bandwidth and a capability for handling aperiodic x^aveforms. 
Kurland & Molony [2] experimentally determined and ranked the error 
performance for 15 different noise waveforms at various signal-to- 
noise ratios (SNR). No pattern appeared to exist in the relative 
standing of the various waveforms as the SNR was changed. Fennick 
[1] has reported on specific switched telephone lines where 2000 dis- 
tinct waveforms have been analyzed. The result was widely varying 
amplitude and phase spectra characteristics. However, when a subset 
of approximately 200 of the waveforms were added together, the com- 
posite noise spectrum appeared to approximate the channel transfer 
characteristic. This would suggest that a weighted unit-impulse 
response of the channel could be used as a typical noise waveform. 

For dedicated lines such as private telephone lines or a Space Shut- 
tle data bus, the channel characteristics will remain constant in 
contrast with the switched telephone system. Furthermore, the number 
of possible noise sources are greatly reduced. Measurements of noise 
which is induced by relay switching have been made on a Space Shuttle 
data bus [33,34] and reveal noise bursts composed of a random number 
of aperiodic waveforms. It is difficult to describe the shape of the 
waveforms from data presented, A reasonable conjecture is that the 
waveforms can be modeled as exponentially-decaying sinusoids. 
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2.2 MATHEMATICAL MODEL 

A typical noise process based on the discussion of Section 2.1 

is modeled as the sum of NB bursts and is shown in Fig. 2-la. The 
till 

k burst, occurs at time and is composed of a random 

number, NG, of aperiodic waveforms as shown in Fig. 2-lb. In a 
restricted special case, a burst consists of only one waveform as 
shown in Fig. 2-lc. 

The generalized impulse noise (GIN) model is formulated as 
NB 

n(t) = ^ b^(t-T^), 0<tlT^, NB = 1,2, (?.l) 

k=l 

It is assumed that each burst, b^(t-T^), k = 1,..., NB, has a finite 
duration TF^, viz., 

b^(t-T^) ^0 \ < t < \+TF^ • (2.2) 


Furthermore, each burst is represented as 


NG 


b. (t) = y 




(2.3) 


j=l 


where NG is the number of sample functions that occur, g, ,(t) is a wave- 

kj 

th 

form that is chosen at random from the k ensemble of sample functions 


(t) , G^(t),..,, G^(t)} , is a random variable weighting function 




, th 


and T, . is the occurrence time of the j sample function. As an exam- 
th 

pie, let the k noise burst contain three sample functions , i . e. , 




c) ExaniDle of Single-Waveform Burst 


Fig. 2-1^ Impulse Noise Model Waveshapes 
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th 

NG = 3, with each chosen from the k ensemble {G^(t), G^Ct), G^(t), 
G^(t), G^(t)} which consists of five sample functions. One possible 
description for the burst is 

Note that the sample functions G^(t), G^Ct), did not appear while 

G^ ( t) was present twice . 

Three specialized versions of the GIN model are now presented. 

Each of these models is formed by restricting the previous model. The 
first specialized version restricts each burst to be composed of sample 
functions selected from the same ensemble {G^(t), G^Ct),..., G^(t)}, 

and will be identified as the GIN (Single Ensemble) model. This allows 
a noise burst to be represented by 

NG 

, 1=1 

where the k subscript for Sj used in (2.3) can be omitted since each 
burst is formed from the same ensemble. The second specialized model 
is obtained by adding the additional restriction that each burst con- 
tains exactly one sample function. Thus, 

NB 

n(t) = ^ I^g^(t-T^) . (2.6) 

k=l 

Note that g^(t-T^^) is randomly selected from the same ensemble of sample 
functions for each value that k assumes. This model will be called the 
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GIN (Single Function) model. The third specialized version restricts 
the sample function of (2.6) to be the same waveform f(t) for each occur 
ence, i.e. , 

NB 

n(t) = ^ 
k=l 

This is identified as the GIN (Unique Waveform) model. If f(t) is 
chosen to represent the channel unit-impulse response then (2.7) is 
similar to the model used by Bellow & Esposito [12] and to the model 
suggested by Fennick [1,20], The idealized impulse used by Ziemer [10] 
and Routs & Moore [11] can be obtained by setting f(t)=6(t). 

Some of the important properties of these noise models are pre- 
sented in Section 2.3. It is shown in Chapter 3 that the performance 
analysis for the GIN (Unique Waveform) model can be extended to include 
the other models as well. 


2 . 3 NOISE CHARACTERISTICS 


Several plausible assumptions will be made concerning the random 
variables involved in the GIN model. Also, the power-density spectrum 
and the p.d.f. for the GIN (Unique Waveform) model will be determined 
as functions of the various noise parameters. 

The reasonable hypothesis that each sample functions, > 

is equally likely to occur at any time in the burst time range, 0 to 
TFj^, leads to a uniform distribution for j j i.e., 
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1 

TF, ’ 

k 


0 < a < TF, . 
— — k 


( 2 . 8 ) 


Consequently 5 for statistically independent occurrence times a Poisson 
distribution with the average occurrence rate v can be assumed for fhe 
number of waveforms, NG, in a burst time TF^, i.e,, 


(v TF, -VTF 

P[NG=m] = j— — e 

m. 


(2.9) 


A similar result applies to the number of noise bursts in time T if the 

n 

bursts are equally likely to occur at any time. 

The p.d.f. descriptions for the random variables associated with 

burst type, sample function selection, and the weighting factor are 

established by the particular system being analyzed. However, it will 

be assumed that the weighting factor p.d.f., p_, (i) , is symmetrical 

about a zero mean. This allows both positive and negative versions of 

the noise waveforms g, ,(t). It is assumed that the statistical descrip- 

kj 

tions of the noise parameters are not time dependent. For some systems 
this will be a valid assumption, while for others it is an approximation 
that must be limited to short time periods. 

It is shown in Appendix A that the GIN (Unique Waveform) model 
yields a one-sided pox^^er-density spectrum given by 


S^(f) = 2 Vga^|F(f)|^ , 


( 2 . 10 ) 


where v is the average occurrence rate of the bursts, a is the vari- 

D 

ance of each weighting factor , and F(f) is the Fourier transform 
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of the noise waveform f(t). The result of (2.10) can be extended to 

the GIN (Single Function) model by taking the expected value with 

respect to the waveform type, cf., (A. 11). Note that the power- 

density spectrum of (2.10) depends on the sample function and is a 

function of frequency, i.e., the spectrum is not necessarily white. 

For the special ease where f(t)=6(t), then F(f)=l and white poise 

2 

results since S (f) = 2v, a . Similar results have been given for the 
n b 

idealized impulse noise model by Ziemer [10] and Routs & Moore [11]. 

The performance results obtained for impulse noise can be coit|- 
pared with those for AWG noise, by relating the dependent variables. 

The convention established for AWG noise is to plot the error proba- 
bilit;y as a function of signal-to-noise ratio, defined as the rptio 
of the average signal energy, to the one-sided noise powep-^ 

density spectrum, N^, i.e., 

g 

SNR = . (2,11) 

o 

Previous results [11] have shown that the error performance for addic- 
tive idealized impulse noise depends on an energy-to-noise parameter 
(ENP), viz.. 


ENP 


^AVG^b 

2 

a 


( 2 . 12 ) 


where is the bit time. Thus for the idealized impulse noise. 


ENP = 2 v„T, SNR . 
b b 


(2.13) 
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For colored Gaussian noise the performance is expressed as a ratio of 
signal power to noise power contained in some specified bandwidth. The 
comparison of SNR and ENP for generalized impulse noise would be more 
complicated than the simple relation shown by (2.13) and xj^ill not be 
determined in this work. 

An interesting characteristic which can be obtained for the GIN 

(Unique Waveform) model is the first-order p.d.f. for the poise n(t), 

i.e., p^(a). Following the lead of Rice [35], Middleton [36], and 

Ziemer [10], it is shown in Appendix A that when the number of noise 

bursts in time T is modeled by a Poisson process with average rate 
n 

Vg, then 


, , ,-1/2 ,(0), , , ^^3^2 .(3), X . 

Pj^(a) = ^2 — 3T~ 


^ 2 ,(4), X . 3^2 ,(6), X 

; ip' ■'(x) ^(x) 


4 : 


72 


(2.14) 


where 


X = V E{I®} f f™(t-x) dr, 
m r> I 

‘0 


(2.15) 


and 


, (k) ^ ^ 1 d^ r -x^/ 2 t , — 

^ (x) = ^ [e ]; X = a/VX 


(2.16) 


dx^ 


- 1/2 

The first term in p (a) is of prder v , the second term of order 

n o 

—1 T“ 3/ 2 

V , and the third term enclosed by the brackets is of order 
B B 

Also, the first term assumes the form of a Gaussian distribution, thus 


the remaining terms show how p (pi) approaches Gaussian as v increases. 

n ' B ' 
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If the Poisson process assumption used to establish (2.14) is 

replaced with the assumption that at most one noise burst can occur 

in time T , the p.d.f. becomes 
n 


p„(») 


= P[NB=0] 6(a) + 


P [NB=1] 
T 


00 X 

If 


Pj(i) 


6 (a-'if (t-t) ) dxdi 


(2.X7) 


In receivers where the bit decision is based on sampling the received 
signal and comparing to a voltage threshold, the p.d.f. ’s of (2.14) and 

(2.17) could be used to determine the probability of error. Using 

(2.17) , the probability that the noise exceeds a positive threshold 
is 


f Pj^M da = J' J pj(i) dxdi, 

TT ■ ^ ‘ -r» 


(2.18) 


where the region is defined such that i*f(t“T^) > The per^ 

formance of a correlation receiver is determined in Chapter 3. 
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CHAPTER 3 


CORRELATION RECEIVER PERFORMANCE ANALYSIS 

The expression for the error performance is derived for a digital 
system that uses a correlation receiver. Although the derivation is 
limited to the GIN (Unique Waveform) model, a procedure is presented 
for extending the results to the unrestricted GIN model. The error 
performance is bounded by a function of the average signal energy, the 
signal set correlation coefficient, and the variance and p.d.f. descra^p-- 
tion of the noise weighting factor. The bound dependence on the p.d.f. 
is removed by applying Chebyshev^s inequality. 


3.1 GENERAL OBSERVATIONS 

A method for judging the error performance of a digital system 
is presented in the section. Subsequently, the analysis for the GIN 
model is obtained in terms of the GIN (Unique Waveform) model. The 
effect of the communication channel on the analysis of the correlation 
receiver is discussed. Several assumptions regarding the GIN (Unique 
Waveform) model parameters are summarised for use in the derivation of 
Section 3.2. 
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3.1.1 Perfomance Criterion 

Two performance criteria x^ere considered for use in this work; 
namely, the probability of bit error and the expected number of bits 
in error given a burst has pccurred. The probability of bit error 
is used extensively fpr AWG noise analysis. There is some precedent 
[6,9-12] for using this criterion in impulse noise analysis; how- 
ever, the more definitive criterion of the expected number of bits 
in error per noise burst has beep employed [20,24] and will be used 
in this work. Each method allows an easy comparison of theoretical 
and experimental results, since the difference between the experi- 
mental numerical average and the theoretical statistical probability 
can be made arbitrarily small as the number of observations is in- 
creased. The expected number of bits in error per noise burst cap 
be converted to the probability of bit error by multiplying by the 
probability of a noise burst and dividing by the expected number of 
bits transmitted during the burst. The relative performance of 
various signaling methods will be based on equal average signal 
energy. 

3.1.2 Rational for Unique Waveform Restriction 

The performance analysis for the GIN model i^ reduced to the 
analysis for the GIN (Unique Waveform) model through a series of 
conditional probability steps. The expected number of bit errors 
per noise burst, E{N}, can be represented a,s a summation of terms 
that are conditioned on the type of noise burst, 
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E{N}= ^ E{N|Burst Type j} P{Burst Type j} . (3.1) 

j 

The evaluation procedure for each conditional term, E{N|Burst Type -j } , 
would be identical to that used for the GIN (Single Ensemble) 
model where only one burst type is allowed. The analysis procedure 
can be further simplified by assuming that the waveforms within a 
noise burst do not overlap. Consequently, the total number of bits 
in error per noise burst, N, is equal to the sum of N^, the number 
of bits in error due to each waveform, viz. 


N = N, +N^+...+N + ...+ N 

12 m NG 


(3.2) 


Furthermore, the expected value of N can be written in terms pf the 
number of waveforms per burst, NG, as 


E{N} 




E{N NG=i} 


P[NG=i] . 


Using (3.2) to express E{N NG=i} gives 


E{N} 




1=1 


P[NG=i] 


V E{N } . 

/— » m 


m“l 


(3.3) 


(3.4) 


The evaluation of E{N } is conditioned on the type of samnle function, 

m y r , 5 

Gk(t) , that occurs and the probability of the occurrence, i.e., 


NE 

E{N } = y 

m ^ 

k=l 


E{N G. (t)} P[G. (t)], 
m K K 


(3.5) 
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where NE is the number of sample functions in the ensemble. Assuming 

that the probability P[G^(t)] is the same for each waveforip in the 

burst, it follows that E{N- } = E{N.}, j = 2,3,.... Consequently, it 

1 i ' 

follows from (3.4) that the expected number of bit errors per noise 
burst is 

NE 

E{N} = i P[NG=i] E{Nj^|Gj^(t) } P[G^(t)]. (3.6) 

1=1 k=l 

This can be written as 

NE 

E{N> = E{nG> Y' E(N^|Gj^(t)} P[Gj^(t)], (3.7) 

k=l 

The evaluation of the terms within the summation would be the sapfie as 
for the GIN (Single Function) model, while the evaluation of each 
E{Ni|Gk(t)} would use the same procedure as for the GIN (Unique Wave- 
form) model. Consequently, a procedure for determining the error per- 
formance of a signal set in the presence of noise described by the GIN 
(Unique Waveform) model can be extended to noise described by the GIN 
model. 

The expected number of bits in error per noise burst which contains 
one unique waveform, i.e., the GIN (Unique Waveform) model, will be 
determined by considering each bit that the noise waveform cojuld pos- 
sibly effect. The data bit which is present at the time of occurrence 

of the noise waveform will be identified as B^. The bits are numbered 

ttl 

consecutively through the K -bit, B , which is the last bit that the 

K. 

noise burst can overlap. An expression for the number of bits in error 


due to the waveform can be written as 
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K 



m=l 


(3.?) 


th 

where x =1 if the m bit is in error and x =0 otherwise. The expec- 
m m 

ted number of bits in error per noise burst, E{N^}, is obtained from 
(3.8) by recognizing that the expected value of a sum of random vari- 
ables is equal to the sum of the expected values. Using the notation 

e{x }=P[e ] i.e., the expected value of x is the probability of 

mm m 

th 

error for the m bit, gives 


K 

E{n} = E{N^} = 

m=l 


3.1.3 Channel and Receiver Considerations 

The digital communication system shown in Fig. 3-1 has been 
selected for analysis. The effect of the linear time-invariant 
channel response, h(t), is resolved prior to the error performance 
derivation in Section 3.2. The channel output, r*^(i) can be written 
as the convolution of the input signal, r(t) and the channel unit- 
impulse response, h(t), i.e., 

r°(t) = r(t)*h(t) = s^(t)*h(t) + n(t)*h(t) = s°(t) + n°(t) . (3.10) 

If the noise is AWG, the minimum probability of bit error is obtained 
by using the f iltered-ref erence correlation receiver [37]. No claim 
of optimality is made when impulse noise is added to the channel. The 
effect of the channel on impulse noise is to change the waveforms in 
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the set {Gt ( t) , G_(t),..., G„„(t)} from which n(t) is composed. Al- 
though the channel changes the signal and noise waveforms, the per- 
formance analysis procedure for the receiver remains the same. Conse- 
quently, the channel response, h(t), will be ignored. 

3.1.4 Assumptions Regarding Noise Parameters 

The GIN (Unique Waveform) model involves four parameters, i.e^, 
the number of noise bursts NB, the unique waveshape f(t), the occur- 
rence time of the burst T, and the weighting factor I. Thq following 
assumptions regarding these parameters are used in Sections 3.2 and 
3.3: (1) the weighting factor I is symmetrically distributed and has 

zero mean, (2) the occurrence time of a burst is uniformly distributed 
over the bit time in which it occurs, (3) the waveform f(t) h^s unit 
energy and the p.d.f. for I is determined from a histogram of the 
square root of the energy in the weighted waveforms , and (4) no more 
than one noise burst can occur in a bit interval. 

3.2 DERIVATION OF GENERALIZED ERROR PERFORMANCE EQUATION 

The digital system shown in Fig. 3-1 and the optimum Gaussian 

correlation receiver of Fig. 3-2 are used in the evaluation of P[e 1. 

m 

The first step is to find an expression for error probability in 
terms of the receiver conditional decision statistic, D^ , j=0, 1, which 
represents the decision statistic, D, given that message symbol m. wa? 
selected at the source. The second sfep uses a sequence of conditional 
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Fig. 3-2. Correlation Receiver 




probability operations to establish a form for the conditional error 
probability. The third step is to normalize the error performance 
equation and simplify the results by introducing the time range func- 
tion 


3.2.1 Error Performance in Terms of Receiver Decision Statistic 

th 

The receiver decision statistics for the m bit which the noise 
waveform affects can be written as 

r(t)[Sj^(t) - s^Ct)] dt - C . (3.11) 

l)T^ 

For the Gaussian noise the bias term C is given by 



C = 





■p[mo]' 


(3.12) 


where , j=0,l, represents the energy of one-si|ied 
noise power-density spectrum. For equal message probabilities the bias 
term reduces to 


C = 




(3.13) 


The assumption that P[mQ] = will be used throughout the remaining 

analysis. Restating (3.11) with the substitution 



24 


and using a corirelation coefficient p defii>ed by 


P 



So(t)Si(t) dt , 


(3.15) 


yields after straightforward manipulation the receiver conditional deci-r 
sion statistic , viz.. 




(t)[s^(t) - 


So(t)] 


dt - 


C. 


(3.16) 


The limits of integration have been dropped and ar^ understood to be 
over the bit interval, and 


Cj = C - (-l)J(pE^-E^) . 


(3.17) 


Substituting (3.13) for C it can be shown that 


S = -^1 = 


(l-2p)E^ + Eq 


(3.18) 


Replacing n(t) in (3.16) with the GIN (Unique Waveform) model of (2.7), 
i.e. , 

NB 

n(t) = I^f(t-T^) , (3.19) 

k=l 

gives 

NB 

^ f(t-T^)[s^(t) - SQ(t)] dt - . 

k=l 


(3,20) 
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Considerable simplification results when it is assumed that at most one 
burst can occur during a bit interval. Conditioning on this assump- 
tion gives 


D^|(NB=1) = J f(t-T^)[s^(t) - s^(t)] dt - C^ , 

I (NB^O) = -Cj , 


(3.21) 


It has been shown [11] that (3.21) can be used to form a lower bound 
to error performance when the number of bursts, NB, is described by 
the Poisson distribution. Furthermore, it was shown that the lower 
bound becomes a good approximation as the product of the average occur- 
rence rate and the bit time T, becomes small, i.e., << 1. 

B b B b 

Considering the receiver decision rule from Fig. 3-2, i.e.. 


D ^ 0 assign output to m^^ 

(3.22) 

D < 0 assign output to m^ , 


it follows that the probability of error conditioned on message symbol 
m. can be expressed as 


P[e |m.] = P[(-1)^D^ >0], j = 0,1 

m ' j j j 


(3.23) 


and the total probability of bit error is given by 

1 

P[e ] = X! P[e |m ] P[m ] . 

™ ^=0 ^ J 


(3.24) 
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3.2.2 Conditional Bit Error Probability Formulation 


Two conditional probability operations are performed on (3.23) in 
order to obtain results which depend on the p.d.f. descriptions for the 
weighting factor and occurrence time.* First the expression involving 
is conditioned on whether or not a burst occurs, i.e., 


P[(-l)^D^ > 0] 



> 0|NB=k] P[NB=k] . 


(3.25) 


For k=0, it is obvious that no error cap occur since from (3.21) 


Thus , 


I (NB=0) = -C^ . 


P[(-1)^D^ > 0] = P[(-1)^D^ > 0|nB=1] P[NB=1] 


(3.26) 


(3.27) 


A normalized probability of error is defined to help simplify the nota- 
tion, i.e. , 


P[e |ni.] . . 

" ~P [fB=i] = > 0 |nB=1]. (3.28) 


The expression for the decision statistic given that one burst has 
occurred is obtained from (3.21), namely 


I (NB=1) = I 


I 


f (t-T) [Sj^(t) - 


SQ(t)] 


dt - 


C. , 
1 


(3.29) 


*An alternative method is presented in Appendix B. 



where the 1 subscripts for I and T have been dropped. The integral 
over the bit interval is represented as 
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V F(T) 


I 


f(t-T)[s^(t) - s^Ct)] dt , 


(3.30) 


where V is the magnitude of the peak voltage of s^(t) and/or s^Ct). 
It follows that 


D^l(NB=l) = I V F(T) - Cj (3.31) 

and 

P, [e |m.] = ?[(-!)’’ I V F(T) > (-1)'^ C.] . (3.32) 

1 m 1 1 

Since Cq=-C^, it follows that the term (-1)^ can be replaced by C^, 
hence 


Pl[£ml®j] = I F(T) > Cq/V] . 


(3.33) 


The right side of (3.33) can be expressed as the marginal probability 
form of a mixed mode probability density function, viz., 


P[(-l) I F(T) > C^/V] ^ 


J j p[(-i)^ 

i T 


i F(t) > CQ/v|l=i,T=T]p^(i)p^(x) drdi . 


(3.34) 


Given values for I and T, the conditional probability term within the 
integrals of (3.34) will equal unity if (-1)^ i F(x) > C^/V and zero 
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otherwise. Thus the integrals over i and t can be restricted to the 
region R. where this inequality is true, yielding 


I™ 

1 m 


.1 -// p,(i 


i)p^(x) dxdi 


R. 

J 


(3.35) 


3.2,3 Normalized Bit Error Probability Formulation 

The bit error probability, P[e^], is obtained by invoking the as- 
sumption that I is symmetrically distributed, i.e,, p^ (i)=p^ (~i) . It 
follows from (3, 34) that 


P[I F(T) > Cq/V] = P[I F(T) < -Cq/V] . (3.36) 


Consequently, substituting (3.33) and (3.36) into (3.24) yields 


1 m 


= [p[mQ] + 


1 m 


mQ] 



Pj(i)p^(x) dxdi . 


(3.37) 


Assuming that T is uniformly idstributed over a bit interval T, , then 

b 

Integrating over x where i F(x) > C^/V gives time ranges T (i) and 
T (i) which depend respectively on the negative and positive values 

Kr 


of i. Hence, 
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Pi [e ] 

1 m 


U 


Tgp(i)Pi(i) di 


oo 

■kl 


tT^(-i) + T^(i)]pj(i) di , 


(3.38) 


where the last equality follows by changing variables and using p^(i) 

= p^(-i). Let F represent the maximum value of |f(t)L and set 
max ^ 115 


I . 
mxn 


V F 


max 


(3.39) 


For values of i < I . both T^„(-i) and T„_(i) are equal to zero. Let 

min BN RP 

T (i) represent the sum of the time range functions. Thus, 

K 

OQ 

^ y Tj^(i)pj(i) di . (3.40) 

^ I . 
min 

It is possible to determine T (i) without first determining T (-i) 

K RN 

and T (i) . This requires that the integration over x in (3.37) be 
performed such that i | F(x) | >C q/V. 

It is convenient to transform I into a normalized (unit variance) 
random variable N. Using the change in variables N = l/o yields a 
density function p^(n) which is equal to ap^(an). Consequently, (3.40) 
can be expressed as 


1 m 


00 


T^(an)p^(n) dn. 


(3.41) 
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The dependence of the time range T (crn) on the variance, can be removed 
if a new function Tj^^(n) is determined by integrating over t such that 
n I F( t) I >Cq/ ( aV) . This gives the desired form for the results as 


^1 

± m 


00 

i I 


TRi(n)Pjj(n) dn 


(3.42) 


I . /a 

min 


The analytical evaluation of T (n) can be extremely difficult; 

Ki 

often, numerical methods are necessary. An illustration of a partic- 
ular method is shown in Fig. 3-3, where the function |f(t) | is plotted 


versus x for 0 < x < T, , 
— — b 


First, a set of K values for n, 


..., n }, are chosen and the corresponding set of slice levels {SL , 
K “ 1 

SL^,..., SL^} are calculated from the expression 
Z K 


SL, = 
k aVn, 


k == 1,2 , ... , K 


(3.43) 


Next, as shown in Fig. 3-3, the plot of |f(x) | is compared to each 

slice level and a set of time range values {T (n-), T (n^),..., 

K1 X Ri Z 

TriCu^)} are calculated. The value obtained in Fig. 3-3 for the 

time range corresponding to SL, is 

k 






(3.44) 


This data can be used to evaluate (3.42) numerically. 
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3.3 PERFORMANCE EQUATIONS FOR SPECIFIC SIGNAL SETS 

The performance analysis which culminated in (3.42) used a general 
form for the signal set (s^Ct), s^(t)}. Although many signal sets have 
been used fpr binary communications, it is possible to classify them 
into a few basic catagories. Unfortunately, a standard identification 
procedure does not exist. For binary signal sets having a correlation 
coefficient, p=0 or p=-l, three clasif ications are proposed as follows: 



Equation 3.48 can be used for Orthogonal signaling (p=0) by replacing 

i|>(t) with (t) . It follows that the time range T^^(n) can be 

determined from the condition |F-(T)| > C^/ [ (1-p ) Van] and that I /a 

I U min 

= ^0^^ ^ where is the maximum value of |f^(T)|. 
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It is possible to express Cq/[(1-p)Vo] in terms of the energy~to- 
noise parameter, ENP, defined by (2.12). First, the average signal 
energy, E , for equally likely message symbols is defined as 

r\V VT 


^AVG 




(3.49) 


For the three classifications under study, C^, as defined by (3.18), 
can be expressed in terms of viz., 


C 


0 


(l-p)E 


AVG “ 


(3.50) 


The voltage V can be written as a function of ENP once the basis func- 
tions, i|^(t) and defined. Two functions are considered: 

namely, ^ (t) for a pulse and (t) for a sinusoid. They are defined 
P ® 

as 


^p(t)| =1 0 1 t < T^ (3.51) 


and 


ijj^(t) = sin(27Tf^t) 0 £ t < T^ . (3.52) 


The corresponding peak voltages, V and V , for On-Off and Antipodal 

p s ^ 

signaling can be expressed as functions of E , i.e., 

AV G 


V 

p 


^ "avg 
(l-p)T^ 


(3.53) 
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and 


= 7"" -V. G , 

V (l-p)T, 


Hence, for these two classes of signal sets 


and 



(l-p)a 


V 20^ 


^ (l-p)ENP 


^0 ^ / (l-p)ENP 

(l-p)o Vg V 4 


th 

It follows from (3.42) that the error performance for the m bit 
given for On-Off or Antipodal pulse signaling as 

00 

TpiCn) Pjj(n) dn , 

-p)ENP/2 
^Imax 

where determined such that n|F^(T)| > /(l-p)ENP/2. For 

sinusoidal basis functions and the same two signal sets the error 
performance is given by 

oo 

\l(n) Pfj(n) dn , 

-p)enpA 4 
^Imax 

where T (n) is determined such that n|F (x)| > / (l-p)ENP/4. In 
Ri i 

similar fashion the result for Orthogonal pulse signaling is 

oo 

- It / ■*" • 

y/ENP 

^Imax 


Pits ] 

1 m 


= kr / 


/a 


i m 




/(T 


(3.54) 


(3.55) 


(3.56) 

is 


(3.57) 


(3.58) 


(3.59) 
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where (n) is determined such that n|F-(T)| > /ENP. For Orthogonal 
R1 J- 

sinusoidal signaling the performance is given by 

oo 

■ kr f 

° /ENP/ 2 
Imax 

where T (n) is determined such that n|F (t)| > /ENP/ 2. 

R1 -L 

Two significant distinctions are found in comparison of the im- 
pulse noise error probabilities with their equivalent AWG noise expres- 
sions. Comparing (3.60) to (3.59) or (3.58) to (3.57) it is seen that 
the lower limits of integration differ by a factor of /?, i.e., the 
choice of basis function to represent the On-Off, Antipodal, pr Orthog- 
onal waveform has an effect on the error performance. This is in 
contrast to AWG noise results where the performance does not depend on 
the basis function. Likewise, comparing Orthogonal to On-Off signaling 
for the same basis function and equal average signal energy, e.g., 
(3.60) to (3.58), reveals a similar difference in the lower limit. 

Again this performance discrepancy is in contrast to AWG results. Con- 
versely, from (3.57) or (3.58) a comparison of the performance expres- 
sions for Antipodal and On-Off signaling indicates the same 3 dB im- 
provement in ENP for Antipodal as was realized in AWG noise theory. 
Performance curves for selected noise waveforms and signal sets are 
presented in Chapter 4. 

3.4 DERIVATION OF PERFORMANCE BOUNDS 

There is motivation for bounding the error performance results 
since determining the noise waveshapes can be extremely difficult. 
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and the results must be specialized for each communication system con- 
sidered. In this section, an error performance bounding procedure is 
developed using an upper bound, , on the receiver decision statistic, 
, In particular V*F(T), cf . , (3.30), is bounded by applying the 
Schwarz inequality, i.e.. 


V-F(T) 


= J f(t-T)[s^ 


(t)-s^(t)] dt 


1 dx J dy. 


(3.61) 


Assuming that the p.d.f. for the weighting factor is determined from 
the square root of the energy contained in the weighted waveforms 
yields 



(m-l)T^^ 


f (t-T) dt 


00 


f^(t) 


dt = 1 


(3.62) 


It follows from substitution of (3.62) into (3.61) that 


V.F(T) £ /(l-2p)E^ + Eq . 


(3.63) 


Consequently, using (3.63) in (3.31) it follows that an upper bound 
can be defined as 


= I /(l-2p)E^ + Eq - Cj ^ I (NB=1) . 


(3.64) 



37 


The probability that > 0 is greater than or equal to the prob- 

ability that (-l)^D^ I (NB=1) > 0. Hence, the error performance, as 
given by (3.28), can be bounded by 


|m.] < P[(-1)^B^ > 0] = P[I > /[(l-2p)E, + E^]/4] . (3.65) 

1 m j — 1 u 

For On-Off, Antipodal, and Orthogonal signal sets it follows from 
(3.50), after changing to the normalized random variable N=I/o, that 


Pl[£n,|mj] 1P[N >V(l-p)E^^g/(2a^)] . 


(3.66) 


This result is independent of m^ , thus 


< 

J 

1 m — 


p^(n) dn . 


ra-p)E 


AVG 


(3.67) 


^ 2a^ 

This bound does not depend on the noise waveform, but the p.d.f. for 
the weighting factor must be known. A bound that is independent of 
p^(n) can be obtained for (3.67) by applying Chebyshev’s inequality, 
which for a zero mean, unit variance, symmetrically distributed ran- 
dom variable states 


P[|n| ^ l] = i P[N ^ l] ^ . 

2L 


(3.68) 


Applying (3.68) to (3.67) yields 
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C 

< J 

1 m — 


Ph<") dP 1 -filpJE 


(3.69) 


f(l-p)E 


AVG 


AVG 


V 2a" 

These results are presented in Fig. 3-4 for the five two-sided p.d.f.s 
of Table 3.1. A comparison of the bounds to performance results for 
specific signal sets is presented in Chapter 4. 


TABLE 3.1 

NORMALIZED PROBABILITY DENSITY FUNCTIONS 


Distribution 


Density Function 


Restrictions 


Gaussian 


Rayleigh 

(Two-Sided) 


Lognormal 

(Two-Sided) 


Uniform 


Hyperbolic 

(Two-Sided) 


— ^ exp[-n^/2] 

Tlir 


n| exp[-n ] 


2/?ir 


a n 


[In |n I H- 
2a^ 


exp 

1 

/i2 


/2 [v/(m-2)(m-3) + Tl Ini]™ 


a > 0 


Ini < »/3 


3 < m < 6 










UPPER BOUND ON P,Ie 
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Fig. 3-4. Performance Bounds 



CHAPTER 4 


SPECIFIC PERFORMANCE RESULTS 


Waveforms are chosen for the GIN (Unique Waveform) model and used 
to calculate numerical performance results, illustrate the effect of 
the weighting factor p.d.f., compare performances of various signal 
sets, and study the tightness of the bound curves derived in Chapter 3. 


4.1 GENERALIZED IMPULSE NOISE WAVEFORMS 

The determination of the impulse noise waveforms present in a 
communication system requires extensive measurements and the use of 
special test equipment* Two specific waveforms, a decaying exponen- 
tial and an exponentially-decaying sinusoid, are selected for closer 
examination, based upon the measurements on the simulated Space Shut- 
tle data bus [33] and the unit-impulse-response conjecture proposed 
in Section 2.1. 

4.1.1 Decaying Exponential 

The decaying exponential is the unit-impulse response of a 
channel modeled by a first-order Butterworth transfer function. 

In order to conserve computer time, the waveform is truncated at 
TF seconds. The expression for the truncated, unit-energy wave- 
form is 


40 
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f(t) = A exp(-bt), 0 £ t £TF (4.1) 

where the decay rate depends on the time constant (1/h) and for the 
assumption of unit energy 

^ """^l - exp(-2hTF) * 

Usually the waveform duration (TF) is chosen to equal several time 
constants. Waveforms are presented in Fig. 4-1 for time constant 
and duration values that are used in obtaining the performance 
results of Section 4.2. A unit-energy pulse is shown for compari- 
son. It follows from (4.2) that the noise waveforms have larger 
peak values as the time constant and duration are decreased. 

4.1.2 Exponentially-Decaying Sinusoid 

The second unit-energy waveform selected is the truncated, 
exponentially-decaying sinusoid, i.e., 

f(t) = C sin(27Tf^t) exp(-dt), 0 t ^ TF (4.3) 

It is assumed that the waveform duration is a multiple of the 
sinusoid period, i.e., 


TF 


n 


f 


o 


n - 1,2,. 


(4.4) 


Hence, to preserve the unit-energy assumption it follows that 
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C = 

The values for d.f , and TF can be chosen to model both baseband 
’ o 

channels without dc response and band-pass channels. 

4.2 PERFORMANCE RESULTS FOR DECAYING EXPONENTIAL 

Several sets of performance results are presented and comparisons 
made to illustrate the effect of the weighting factor p.d.f., signal 
set, and noise waveform. Results for the five signal sets shown in 
Fig. 4-2 with a 1 ms bit time are obtained using the Fortran-IV com- 
puter program GINIMP which is described in Appendix C. A frequency 
of 1 kHz is used for the ASK and PSK signal sets. For FSK, a fre- 
quency of 1 kHz is used for s^(t) while 2 kHz is used for s^Ct). The 
program requires approximately 2.5 minutes of IBM System/ 360 computer 
time for each overlapped data bit and ENP value analyzed; consequently, 
only three or four points are computed for each case and a smooth 
curve drawn through the points. to obtain the performance curve. Bound 
curves are displayed for comparison with the performance results. 

These bounds are obtained after multiplying the value obtained 

from Fig. 3-4 by the number of data bits that the noise waveform can 
possibly overlap. Although the bounds are not tight, they do have the 
general shape of the actual performance curves and can be obtained with- 
out knowledge of the noise waveform. 
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4.2.1 Weighting Factor Distribution 

The effect of changing the weighting factor (I) p.d.f. is 
demoTist rated in Fig. 4-3. Unipolar signaling and the decaying- 
exponential noise waveform with a 0.2 ms time constant are used 
to obtain performance results for the five weighting factor p.d.f.s 
of Table 3.1. The Gaussian, Rayleigh, and uniform distributions 
yield results which are approximately equal. The results for the 
lognormal and hyperbolic distributions were obtained for parameter 
values a=0.1 and m=3.5, respectively. The parameter value controls 
the shape of these p.d.f.s and should have a corresponding control 
in the shape of the error performance curves. Results for typical 
parameter values in these distributions and an idealized impulse 
waveform have been presented by Routs St Moore [11]. The Gaussian 
p.d.f is selected as being typical and is used exclusively in 
determining the remaining performance results. 

4.2.2 Waveform Time Constant 

The effect of changing the decaying-exponential time constant 
is studied in Fig. 4-4 for Unipolar signaling and a Gaussian p.d.f. 
It was necessary to increase the waveform duration (TF) for the 
1 ms time constant in order to calculate its true impact on the 
data. The expected number of bit errors increases with the time 
constant, which is intuitively satisfying since for large time 
constants the noise waveform more closely approximates the Unipolar 
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1 



-50 -40 -30 -20 


ENP (dB) 


Fig. 4-3. Effect of Weighting Factor P.D.F. on Unipolar 

Signaling in the Presence of Decaying-Exponential 
Noise with 0.2 ms Time Constant 
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Fig. 4-4. Effect of Decaying-Exponential Waveform Time 
Constant on Performance of Unipolar Signaling 
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waveform. Hence, the correlation receiver would be more likely to 
mistake the noise for a signal bit. 

A similar presentation is made in Fig. 4-5 for ASK signaling. 
Note that for ENP values greater than -43 dB, the curve for the 1 
ms time constant indicates a lower error rate than the curves for 
the shorter time constants. A plausible explanation for this 
apparent contradiction to the trend is the lower peak cross- 
correlation which results from the unit-energy noise waveforms. 

The peak value of the noise must decrease as the time constant 
and/or duration are increased. As the peak value decreases, the 
minimum weighting factor, I . , must increase before an error 
occurs; however, the probability of I assuming values greater 

than the new I . is smaller, thus the error rate decreases. The 
min ’ 

effect of the peak values in the correlation function becomes evi- 
dent as the ENP value increases. The bound curves presented in 
Figs. 4-4 and 4-5 are obtained by doubling (3.69); hence, they 
only apply to the noise waveforms of 1 ms duration. The bound 
for the 4 ms waveform appears in Fig. 4-10. 

The* relative performance of the Unipolar and ASK signal sets 
is dependent on the time constant of the decaying-exponential noise 
waveform, as illustrated in Fig. 4-6. The ASK performance shows a 
5 to 10 dB improvement for the larger time constant. This reflects 
the generally lower ASK-noise crosscorrelation values (including 
peak values). For the 0.1 ms time constant, the Unipolar-noise 
crosscorrelation function (including peaks) is slightly larger than 
for ASK signaling, which explains the better performance of ASK at 
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Fig. 4-5. Performance of ASK Signaling for Various 
Decaying-Exponential Time Constants 
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Fig. 4-6. Comparison of Unipolar and ASK Signal 
Sets for Decaying-Exponential Waveform 


51 


low ENP values. However, the ASK signal set has a larger peak value 
than the Unipolar set, which is reflected in the performance equations 
of (3.57) and (3.58), and accounts for the crossover of the curves 
for larger ENP values. 

4.2.3 Signal Set 

An error performance comparison for five signal sets is shown 
in Fig. 4-7 for a Gaussian p.d.f. description of I and the decaying- 
exponential noise waveform with a 0.2 ms time constant. Several 
conclusions can be made from a study of these results. The error 
performance of Polar (or PSK) signaling exhibits a 3 dB improvement 
in ENP value over that of Unipolar (or ASK) . This result agrees 
with that observed for AWG noise. Conversely, it is apparent that 
a difference exists in the performances of Unipolar and ASK signals 
(or Polar and PSK) which is not present with AWG noise. Furthermore, 
the FSK signal set performance appears to approach that of ASK for 
large ENP values, which contradicts the AWG noise results for which 
the signal sets have equal performance. Furthermore, for ENP values 
below -33 dB the FSK error rate becomes less than that for PSK. 
Although the FSK-noise crosscorrelation is generally less than that 
obtained with the single sinusoid used for ASK and PSK, it does have 
slightly larger peak values. This explains the change in the relative 
performance of FSK as the ENP value increases. 

The bound curves in Fig. 4-7 are upper bounds and it appears that 
the bound for p = -1, i.e.. Polar or PSK signaling, could be used for 
all five signal sets. However, this conclusion is not valid, as it 
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Fig. 4-7. Comparison of Five Signal Sets for Decaying- 

Exponential Waveform with 0.2 ms Time Constant 
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applies only for this particular set of assumptions regarding noise 
waveform, time constant, and weighting factor p.d.f. 

4.3 PERFORMANCE RESULTS FOR EXPONENTIALLY-DECAYING SINUSOID 

The effect of changing the noise waveform is demonstrated by 
considering the exponentially-decaying sinusoid of (4.3). Unipolar 
and ASK signal sets are compared in Fig. 4-8 for two values of the 
noise time constant and a Gaussian p.d.f. The higher error rate for 
ASK is explained by the larger correlation between the signal and 
noise waveforms. Similarly, the improvement experienced by Unipolar 
signaling as the time constant is increased is due to the decreased 
correlation between signal and noise. The relative performance of 
the two ASK curves reflects a generally larger crosscorrelation with 
the noise for the 2 ms time constant, but a larger peak cross- 
correlation for the 1 ms time constant. 

A comparison of the five signal sets is presented in Fig. 4-9 
for the exponentially-decaying sinusoidal noise waveform with a time 
constant of 1 ms. These results indicate the same 3 dB improvement 
of Polar over Unipolar and PSK over ASK as is obtained with an AWG 
noise assumption. The FSK performance is almost the same as the PSK 
over the range of ENP values shown, and thereby differs from the 3 dB 
differential for AWG noise results. The performance difference 
between Unipolar and ASK also contrasts with AWG results. 
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Fig. 4-8. Comparison of Unipolar and ASK Signal Sets 
for Exponentially-Decaying Sinusoid 
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Fig. 4-9. Comparison of Five Signal Sets for Exponentially- 
Decaying Sinusoid with 1 ms Time Constant 
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4.4 EFFECT OF SIGNAL AND NOISE CORRELATION 

The performances of Unipolar and ASK signals are compared in 
Fig. 4-10 using both the decaying-exponential and the exponentially- 
decaying-sinusoidal noise waveforms with the same time constant and 
a Gaussian p.d.f. for I. It is noted that the performance of ASK 
with the decaying-exponential noise is approximately equal to that 
of Unipolar with exponentially-decaying-sinusoidal noise. This is 
reasonable, since in both cases the signal and noise crosscorrelation 
are essentially equal. The larger peak value of the ASK signal set 
dictates the crossover in the relative performances as the ENP value 
increases. When the noise waveforms are interchanged, the ASK sig- 
nal set has a small advantage over Unipolar, although both perfor- 
mances are inferior to the aforementioned combinations. 

4.5 COMPARISON TO RESULTS FROM THE LITERATURE 

The relative performance of various signal sets has received 
some attention in the literature. Usually, the comparisons are 
between signaling methods that utilize a carrier, e.g., ASK, PSK, 
or FSK. For idealized impulse noise and coherent demodulation, 
it has been shown by Ziemer [10] and Houts & Moore [11] that PSK 
exhibits a 3 dB improvement over ASK. Houts and Moore have also 
shown that the performance of FSK differs from ASK. Engels [24] 
and Bennett & Davey [38] have used more realistic noise models to 
show that PSK performs better than FSK. Experimental measurements 
reported by Bodonyi [23] have shown the presence of crossovers in 
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Fig. 4-10. Comparison of Unipolar and ASK Signal Sets 
for Waveforms with 1 ms Time Constants 
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the performance curves for noncoherent FSK and coherent PSK sig- 
naling techniques. However, for large SNR values it appears that 
PSK has an advantage. These results for carrier systems appear 
to support Bennett & Davey’s hypothesis that the relative error 
rates of data transmission systems can be determined for nearly 
all noise environments by comparing performances in the presence 
of AWG noise. 

Most of the above results are consistent with the results 
of Sections 4.2, 4.3 and 4.4. However, the performances of the 
On-Off, Binary Antipodal and Orthogonal signaling methods were 
shown to be dependent on the signal-set basis function, e.g., 
pulse or sinusoid. This represents a major exception to the 
AWG-relative-performance hypothesis proposed by Bennett & Davey. 
This discrepancy is caused by the error performance for impulse 
noise being dependent upon the signal-to-noise correlation, 
whereas AWG noise is uncorrelated with all signal sets. New 
conclusions have also been obtained. The 3 dB improvement of 
PSK over ASK has been established for the GIN (Unique Waveform) 
model, and for the conditions investigated, the FSK signal set 
can perform better than PSK for low ENP values. Similarly, a 
study of pulse signaling has revealed that the 3 dB advantage 
exist for Polar over Unipolar. 



CHAPTER 5 


METHODS FOR IMPROVING RECEIVER PERFORMANCE 

The linear correlation receiver does not represent the optimum 
receiver for impulse noise corrupted signals. In fact, a variety of 
methods have been proposed for improving the receiver design. A sur- 
vey of these recommendations is presented in Section 5.1, after which 
a specific improvement method is selected for analysis and the perfor- 
mance results are compared to the linear correlation receiver. Many 
factors other than the receiver design can affect the performance of 
a digital system, e.g., noise source proximity, system shielding, 
ground loops, etc. Such factors are not considered here. 

5.1 SURVEY OF RECOMMENDATIONS FOR IMPROVED RECEIVER DESIGNS 

The topic of designing receivers for impulse-noise corrupted 
signals has not received extensive treatment in the literature. 
However, several intuitively-motivated techniques have been reported 
which employ theoretical and/or experimental procedures. Generally, 
the proposed improvements attempt to counteract either the large 
amplitude or aperiodic nature of the impulse noise and will be clas- 
sified as either amplitude- or time-based modifications to existing 
receivers . 

The amplitude-based modifications act to suppress the large 
amplitude excursions of the impulse noise by including a nonlinear 
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device in the receiver. As discussed in Section 1.2, Snyder [31] and 
Rappaport & Kurz [32] have used the likelihood ratio to investigate opti- 
mum receivers for specific impulse noise models. Each design resulted 
in the use of a nonlinear device to suppress large amplitude excursions 
in conjunction with a linear receiver. The specification of the non^ 
linear device depends on the noise statistics and must be changed if 
the noise description changes. Sub-optimum receiver designs using 
nonlinear devices have been proposed with either intuitive or empiri- 
cal justification. Specifically, several authors [13,26,27,39,40] 
have proposed clipping devices which limit the maximum magnitude of 
the received signal. Another suggestion [29] is to use a blanking 
device which has a zero output when the input signal exceeds specified 
positive and negative thresholds. It has been postulated that the 
reduction in noise power achieved through the use of clipping or blank- 
ing devices yields improved error performance; however, there is no 
such guarantee because of the possibility of improving the signal-to- 
noise cross correlation. Bello & Esposito [12] have reported the per- 
formance improvement obtained using a clipping device with PSK signal- 
ing and an atmospheric noise nodel. They assumed that the amplitude 
of the noise was much greater than that of the data signal. 

The time-based modifications exploit the aperiodic and frequency 
spectrum properties of the impulse noise. Examples of this technique 
are the smearing-desmearing [30] and swept-f requency-modulation [41,42] 
systems. These systems use a special modulator or a filter to destroy 
the phase relationship of the components in the transmitted signal. At 
the receiver the inverse of this operation is performed which yields 
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the original phase relationship for the data signal, but destroys the 
phase relationship of the channel noise components* These techniques 
have not been employed extensively because of their poor performance [43], 
One possible explanation for the poor performance is that the true char- 
acteristics of impulse noise are not represented by a model where ampli- 
tude peaks are caused by a phase relationship of noise components. 
Another time-based modification uses a pulse-width discriminator [28] 
to decide if the received signals are noise or data pulses. Those that 
qualify as data bits produce fixed-width pulses at the output of the 
discriminator. This system can only be effective if the time duration 
of the noise waveforms are short, compared to the data bit. 

5.2 PERFORMANCE ANALYSIS FOR A NONLINEAR CORRELATION RECEIVER 

The error performance of the system shown in Fig. 5-1 is deter- 
mined using the general approach developed in Chapter 3. The first 
step is to find the conditional decision statistic, , for the m^^ 
bit, viz.. 


where { represents the operation performed by the nonlinear device. 
Equation 5.1 can be changed into the form of (3.16). First, is 

added and subtracted to the { term and the following notation is 

introduced. 


mT- 

/ 

(m-l)T. 


(t) + n(t)}° [s^(t) - SQ(t)] dt - C , 


(5.1) 


{n(t)}°° = {Sj (t) + n(t)}° - s j (t) . 


(5.2) 
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Fig. 5-1. Digital System with Nonlinear Correlation Receiver 
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Incorporating (5.2) into (5.1) and dropping the limits of integration 
gives 


= J* {n(t)}°° 


[s^(t) - SQ(t)] dt - , 


(5.3) 


where is defined by (3.17). Assuming noise which is represented by 
the GIN (Unique Waveform) model and ignoring the possibility of more 
than one burst in a bit interval yields 


D^|(NB=l) = J {1 f(t-T)}°° 


[Sj^(t) - 


SQ(t)] 


dt - 


C. 

3 


(5.4) 


The normalized random variable N=I/a can be introduced by defining 
another nonlinear operation, namely, 


a {N f(t-T)}°°° = {a(l/o) f(t-T)}°° . 


(5.5) 


Thus 


D^|(NB=1) = a / {N f(t-T)} 


/ 


lOOO 


[s^(t) - 


SQ(t)] 


dt - C. 


(5.6) 


Factoring the peak voltage V from the signal set and defining 


V 


Fj (N,T) 


= / 


{N f(t-T)} 


OOO 


[s^(t) - 


SQ(t)] 


dt 


(5.7) 


yields 


I (NB=1) = a V Fj (N,T) 


C. . 
3 


(5.8) 
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Following the procedure used in obtaining (3.32) the normalized proba- 
bility of bit error can be written using (5.8) as 


Plt£ml“j3 = P[(-l)^ F^(N,T) > (-1)^ C^/(aV)] 


(5.9) 


Alternatively, (5.9) can be represented as 


PlESmlm.] = 


= JJ P[(-l)^ Fj(n,x) > (-1)^ Cj/(aV) |N=n,T=t] 


n T 


Pjj(n)p^(T) dTdn , 


(5.10) 


or 


Plt^^lmj] 


• //-. 


(n)p^(x) dxdn , 


(5.11) 


where designates the region of N and T such that the inequality 
established on F^(n,x) in (5.10) holds. Assuming that T is uniformly 
distributed over the bit interval and integrating over t yields 

-N 

PlEEjnlnij] J" T^(n)pj^(n) dn + ^ J~ T^(n)p^(n) dn , (5.12) 

^ ~0O ^ N 


where the time range T^(n) equals zero for -N^<n<N^. 

The nonlinear device must be specified before numerical results can 

be obtained. The limiter model of Fig. 5-2 represents an extension of 

the clipping and blanking concepts discussed in Section 5.1 and is used 

by the computer program LIMITER, which is described in Appendix D, to 

determine the error performance given by (5.12). Values for a and a 

p n 




Fig. 5-2. Limiter Model 


correspond to the magnitudes of the peak positive and negative signal 
set voltages, respectively. For ASK, PSK, FSK and Polar signal sets 


a = a , while for Unipolar and other one-polarity signal sets, a =0. 
p n ° n 

The values of 3 and 3 are parameters specified by the user, thus 
P n 

either a specific limiter can be modeled or the values changed to 
determine their effect on the error performance. 


5.3 COMPARISON OF LINEAR AND NONLINEAR CORRELATION RECEIVERS 


Performance results for the linear and nonlinear correlation 
receivers are compared for several signal sets. The results for the 
nonlinear device were obtained by running the LIMITER program on the 
UNIVAC-1108 computer. Each value calculated required approximately 
10 minutes of CPU time for each data bit analyzed, i.e., 40 minutes 
of computer time was required to obtain one point if the noise wave- 
form overlapped four data bits. For this reason, a comprehensive 
study of signal sets, weighting factor p.d.f. and noise waveform param- 
eters was not attempted. The results presented assumed a Gaussian 
p.d.f. for the weighting factor. 

The error performance curves presented in Fig. 5-3 illustrate the 

effect of changing the clip levels 3„ and 3 . When 3 and 3 =0, a 

np PPn’ 

performance improvement of approximately 6 dB is obtained in comparison 
with the case where This is intuitively satisfying, since 

setting 3^=0 is a logical choice for Unipolar signaling. For both sets 
of clip-level values the performance improvement over the linear corre- 
lation receiver decreases as the ENP value increases. However, at a 
-3 

10 error rate there is still an impressive 5 dB improvement over the 
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Fig. 5-3. Effect of Nonlinear Receiver on Unipolar Signal 

Performance in the Presence of Decaying-Exponential 
Noise with a 0.2 ms Time Constant 
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linear receiver for the case when 3 =ot and 3 =0, It is reasonable 

P P n 

to assume that the performance improvement decreases as ENP increases 
since the data signal begins to dominate the noise and the clipping 
device will have less overall effect. A similar conclusion is obtained 
from a study of Fig. 5-4 for ASK signaling and the decaying-exponential 
noise waveform. The clip levels correspond to the peak voltage of the 
signal set, i.e. , before, the performance improvement 

decreases as the ENP value increases. It was shown in Chapter 4 that 
the Unipolar signaling performance was poorer than ASK when linear cor- 
relation detection was used in the presence of decaying-exponential 
noise. This disparity was attributed to the higher crosscorrelation 
between the Unipolar signal and the noise. However, a comparison of 
the error performances shown in Figs. 5-3 and 5-4 indicates that the 
Unipolar signal performs better than ASK when the limiter is inserted, 
at least for the chosen values of 3 and 3 . 


p n 

The exponentially-decaying sinusoidal noise waveform with a 2 ms 

time constant was used to investigate the relative performance of ASK 

and PSK signaling. A comparison of ASK signaling performances for the 

linear and nonlinear receivers is shown in Fig. 5-5. The clip levels 

were not chosen to equal the peak signal value, instead they correspond 

to the rms value, i.e., 3 =3 =a lyfl. Although for low ENP values a 

very small improvement is obtained, the curves have converged for an 

ENP greater than -30 dB. The PSK signaling performance is shown in 

Fig. 5-6 for two versions of the limiter. When 3 =3 =a /vT, the per- 

P ^ V 

formance is actually poorer than it was for a linear correlation detec- 
tor. It is evident that improved performance can be obtained through a 
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Fig. 5-4. Effect of Nonlinear Receiver on ASK Signal Performance 
in the Presence of Decaying-Exponential Noise with a 
0.2 ms Time Constant 
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Fig. 5-5. Effect of Nonlinear Receiver on ASK Signal Performance 
in the Presence of Exponentially-Decaylng-Sinusoidal 
Noise with a 2 ms Time Constant 
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Fig. 5-6. Effect of Nonlinear Receiver on PSK Signal Performance 
In the Presence of Exponentlally-Decaylng-Slnusoldal 
Noise with a 2 ms Time Constant 
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more logical choice of clip levels, i.e., however, the effect 
is minimal for either set of clip levels. Consequently, with a limiter, 
PSK does not possess the full 3 dB advantage over ASK which was present 


for linear correlation detection. 


The results presented in this section indicate that error perfor- 
mance can be improved by the proper choice of clip levels in the non- 
linear receiver. The improvement obtained by using the limiter would 
probably become more pronounced as the noise time constant decreases, 
since the peak noise voltage increases to meet the unit-energy assump- 
tion for the basic waveform f(t). 



CHAPTER 6 


CONCLUSIONS AND RECOMMENDATIONS 


6.1 CONCLUSIONS 

The generalized impulse noise (GIN) model can approximate the 
noise burst structure observed in many digital systems provided the 
average burst rate, expected number of waveforms per burst, sample 
function ensembles associated with the burst, and the weighting 
factor p.d.f. can be obtained. A specialized model, the GIN (Unique 
Waveform) model, is readily obtained by assuming that each noise burst 
contains a randomly weighted version of the same basic waveform. It 
has been shown that the power-density spectrum for this specialized 
model is not white because of the dependence on the noise waveform. 

The GIN (Unique Waveform) model is used to determine the error 
performance equations for a correlation receiver. These results 
express the expected number of bit errors per noise burst as a func- 
tion of the following parameters: (1) energy-to-noise parameter (ENP) 
of (2.12), (2) noise amplitude-weighting-factor p.d.f., (3) signal- 
set correlation coefficient, and (4) a time range function which 
depends on cross correlation between the signal-set basis functions 
and the noise waveform. The performance analysis for the GIN (Unique 
Waveform) model can be extended to the general model. 

Comparison of the performance results for impulse and additive 
white Gaussian (AWG) noise reveals two major differences. For 


73 



74 


impulse noise, Orthogonal signaling does not yield the same error 
performance as On-Off signaling with the same average energy. Also, 
the choice of the signal-set basis functions has an effect on the 
performance. These conclusions are supported by the dependence of 
impulse noise performance results on the correlation between the 
signal and noise waveforms. A conclusion which agrees with the AWG 
result is that Antipodal signaling has a 3 dB advantage over On-Off 
signaling. 

The error performance can be upper bounded by equations which 
depend on the average signal energy, signal-set correlation coef- 
ficient, plus the noise weighting-factor variance and p.d.f. Spe- 
cifically, the bounds do not depend on the signal-set basis functions 
or the noise waveform. Another upper bound, which is also independent 
of the weighting factor p.d.f., is formed by applying Chebyshev’s 
inequality to the previous bound. It has been shown that for impulse 
noise the linear correlation receiver is not optimum, i.e., the proba- 
bility of error is not minimized. This conclusion is supported by the 
dependence of the performance equations on the signal and noise cross- 
correlation. One method for improving the error performance is to 
insert an amplitude-dependent limiter, i.e., a nonlinear device, prior 
to the receiver input. It was shown that the error performance is 
dependent upon the choice of parameters used to specify the limiter. 


6.2 RECOMMENDATIONS FOR FURTHER STUDY 

Further work can be performed in several areas of this inves- 
tigation. The relationship between signal- to-noise ratio (SNR) and 
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ENP has been established only for white Gaussian noise and idealized 
impulse noise. It would aid the comparison of results if the param- 
eters were related for a bandlimited channel and the GIN model. This 
study has concentrated on the analysis of the correlation receiver, 
but an investigation of sampling receivers could be performed as sug- 
gested by the procedure of Chapter 2. The extension of the nonlinear 
receiver designs of Chapter 5 to an optimum receiver for impulse noise 
offers an opportunity for a significant contribution. 

Although the bounds of Chapter 3 are useful, other techniques 
should be investigated in an attempt to improve the tightness of the 
bounds. It might be possible to approximate the error performance of 
the GIN (Single Ensemble) model by using an average waveform G(t) which 
is obtained by averaging over the sample functions in the ensemble. 

For example, the error expression of (3.7) might become, 

NE 

E{N} = E{NG} E{N^|Gj^(t)} P[Gj^(t)] E{NG} E{Nj^|G(t)} (6.1) 

k=l 


where 


NE 

G(t) = Y P[G^(t)] (6.2) 

k=l 

Once G(t) is obtained, the evaluation of E{N^|G(t)} is equivalent to 
the techniques developed herein for the unique waveform model. Finally, 
the derivation of upper bounds for various nonlinear correlation receiv- 
ers would reduce the computer time required to obtain specific results. 



APPENDIX A 


DERIVATION OF POWER-DENSITY SPECTRUM AND PROBABILITY DENSITY 
FUNCTIONS FOR THE GIN (UNIQUE WAVEFORM) MODEL 

In order to characterize the GIN (Unique Waveform) model it is 
helpful to specify both the noise power-density spectrum and ampli- 
tude probability density function (p.d.f,)* These expressions are 
derived in the following sections. In particular, two p.d.f. descrip- 
tions are obtained for the noise model. The first p.d.f. is based on 

the assumption that no more than one noise burst occurs in time T , 

n 

and the second results from using a Poisson distribution to model the 
number of noise bursts. 


A. 1 POWER-DENSITY SPECTRUM 


The power-density spectrum can be determined from the expression 
for the average power, by 


P 

AVG 



(A.l) 


where T^ is the duration time of the noise and the expected value E{ } 
is taken with respect to all of the random variables involved in describ- 
ing the noise process. The first step is to apply Parseval’s theorem, 
vi?. , 
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1 

T 

n 



n^(t) dt 


1 _ 

T 

n 


oo 

J 

— oo 


df 


(A. 2) 


where N(f) is the Fourier transform of a sample function from the ngise 
process. For the GIN (Unique Waveform) model of (2.7), 

NB 

N(f) = ^ ^ (A. 3) 

k=l 


where F(f) is the Fourier transform of f(t). Thus 

NB NB 

|N(f)|^ = N(f).N*(f) = ^ Ij^Ij^|F(f) 1^ . (A. 4) 

k=l i=l 

Substituting (A. 2) into (A.l) and interchanging the order of integra- 
tion and expectation gives 


AVG 


oo 

= ^ f E{|N(f) 
n 


■} df . 


(A. 5) 


2 

The expected value of |N(f)| with respect to I is given by 


NB NB 


Ej{|N(f) 


} = ^ e j^TTfdk Ti) ^ 


k=l i=l 


For the case where and are zero-mean independent random variables , 

2 2 
it follows that ^ where 6^^ = 0 if i ^ k and is 

2 2 

the weighting factor variance. If is a constant a for all bursts, 
then (A. 6) becomes 


Ej{|n(F)|^} = NB |F(f)|^ 


(A. 7) 
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Taking the expected value of (A. 7) with respect to NB yields 


E{lN(f)|2} = Ej^g{E^{|N(f)|^}} = 0^ |F(f)l2 E^g{NB} . (A. 8) 


For a process where E^^{NB} = (e.g., a Poisson process with 

average rate v ) the average power expression of (A. 5) becomes, 

B 


AVG 


oo 

■ / V 


^ 'F(f) 1^ df 


(A.9) 


It follows from (A.9) tha^ the one-sided power-density spectrum is 


S^(f) = 2 VgO^ |F(f)|^ 


(A. 10) 


This result can be extended to the GIN (Single Function) model of 
(2.6) by replacing f(t) with the sample function G^(i) taking 

the expected value with respect to the type of waveform. The result 
is 

NE 

S^(f) = 2 VgO^ ^ P[G^(t)] , (A.ll) 

k=l 

where ^^^{G^(t)} is the Fourier transform of G^(t). 


A. 2 PROBABILITY DENSITY FUNCTIONS 

The GIN (Unique Waveform) model is used in determining p^(ot), 
the amplitude p.d.f. for impulse noise. Two random variable models 
are examined for NB, the number of bursts that occur in time T . 
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The procedure used in the derivation is similar to that employed by 
Rice [35], Middleton [36], and Ziemer [10], The characteristic func- 
tion M(u) is determined first and the noise p.d.f. is obtained through 
an integral transform. 


A. 2.1 Basic Form of Characteristic Function 


The characteristic function M(u) can be expressed in temns of 
the number of noise bursts that occur, viz.. 


M(u) = M(u|nB) P[NB] , (A. 12) 

NB=0 

where 

NB 

M(u|nB) = |nb} = E^exp[ju ^ l^f(t-Tj^)]| . (A.13) 

k=l 

It is assumed that the random variables and T^ are statistically 
independent and that L (or T, ) is independent of I. (or T.), j ^ k. 
It is further assumed that the p.d.f. descriptions for the weighting 
factor and occurrence time do not change between bursts. Hence, the 
subscripts are dropped and 


M(u|nB) = T)j^NB 


(A. 14) 


The expected value is given by 
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T)j ^ J JejuifCt T) 
i T 

Thus, the characteristic function can be expressed as 


M(u) 


oo 

■ 


julf (t-T) 


tNB 


Pj(i) 


P^(t) 


dxdi 


P[NB] 


(A. 16) 


Two models for P[NB] are examined, namely a single-burst assumption 
and multiple bursts with occurrence times described by a Poisson dis- 
tribution. 


A. 2, 2 Single-Burst Assumption 

For some systems it is appropriate to assume that no more than 
one noise burst can occur in the time range of interest, i.e,, 

P[NB ^ 2] = 0. With this assumption (A, 16) becomes 


M(u) = P[NB=0] + P[NB=1] J' J P^-Ci) P^(t) dTdi . 

i T 


(A, 17) 


The corresponding p.d.f. can be obtained from the integral transform 


P„(o) 


oo 

.i- r 

2tt J 


M(u) 


e ^ du 


(A. 18) 


Substituting (A. 17) into (A. 18) and interchanging the order of inte- 
gration gives , 


■ n 


p^(ot) = P[NB=0] 6(a) + P[NB=1]| / 6 (a-if ( t-x) )p^ (i)p^ (x) dxdi. (A, 19) 


X X 
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For given Pj(i) this expression is difficult to evaluate 

because of the delta function 6(a-if (t-x) ) , but it can be used readily 
in determining functions that involve the integral pf p^(a), e-g., the 
cumulative distribution function. 


A. 2. 3 Poisson Distribution Assumption 

The characteristic function of (A. 12) can be expressed using a 

Poisson distribution to describe the number of noise bursts in time 

T , i.e. , 
n’ * 


M(u) = ) M(u|nB) 

NB=0 


NB 


nb: 




(A. 20) 


An expression for M(u|nB) can be obtained by expanding the exponen- 
tial term in (A. 14) into a power series and taking the expectation 
of each term, namely, 


m(u|nb) 


■ [e 

m=0 


(ju) 


T NB 


E{I‘“} E{f“‘(t-T)} 


(A. 21) 


where 


E{l“} 


J~ i“ Pj(i) di 


(A. 22) 


and 


E{f"^(t-T)} 


OO 

J" f^Ct-x) p^(t) dx 

— OO 


(A. 23) 



82 


If the burst occurrence times are uniforaly distributed over T^, then 


E{f“(t-T)} 


n 


f^(t-T) dT 


(A. 24) 


Substituting (A. 24) and (A. 21) into (A. 20), gives 


M(u) 


TvTTj— o L n 


E{I‘“} J f”"(t-T) dT 
0 


NB. I a m; 

NB =0 *- m =0 


NB 


. (A. 25) 


Letting 


X = E{l"*} f f“(t-T) dT , 

in D I 

0 


(A. 26) 


then 


M(u) = e 


oo oo 

■ I* [I 

NB =0 ^ m =0 


(iu) 


m 1 NB 


m m* 


= e 


-Vt.T. 


B^n 


exp 


oo 

ll 

m==0 


(ju) 


m 


m m: 


exp 


oo 

[I 

m=l 


(A- 27) 


When I is symmetrically distributed about a zero mean, it follows that 
E{i} = 0 and from (A. 26) that = 0. The p.d.f. for the noise can be 
determined from the transform relation of (A. 18). It follows that 


oo 

- h J 


exp[-jua-u X^/2] exp 


[^■^1 


du 


(A. 28) 



A series solution can be obtained by a two-step procedure, namely 
replace the second exponential by 
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7 X + ... . (A.29) 

I Z—j m ml I Z— I mm. 12 ^ 

m=3 m=3 


exp 


and perform a term-by-term integration using the integral relation 


oo 

r exp[-jua-u^X /2] du = (-1)” - exp[-a^/2X„] . (A. 30) 

J ^ ^ h da^ ^ 


The derivation of (A. 30) is presented in Section A. 2. 4. The p.d.f. 
becomes 




-2 


. N 1 .(0)/ , 3 2 .(3), , 

P„(“) = — M — <P (x) 


" ^ 


-5/2 


XX ''''' X^X 

. .42 ,(4), , . 3'^2 ,(6), , 

+ I TJ ^ (x) + ^2 


41 


+ • 


(A. 31) 


where x = a/t/^ and 


»(n)„ 

/?ir dx 


-x^/2 


(A. 32) 


Note that terms have been collected according to the combined powers 
of the A terms, or equivalently to the power of v^, cf., (A. 26). Since 
the first term is of the Gaussian form, it is evident that the p.d.f. 
approaches Gaussian as becomes large. 
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A. 2.4 Derivation of (A. 30) 

The integral relation of (A. 30) is obtained from the definite 
integral 


r -V- , 

/ e - dy 


= /rr 


(A. 33) 


Changing variables y = t v4i/2 yields 


oo 

J dt = /2Tr/h 


(A. 34) 


Differentiating n times with respect to h and multiplying by (-1)^ 2^ 
gives 


oo 

J t'" dt 

_ OO 




(2n)I ,-n-l/2 
n 

Z n. 


(A. 35) 


Next, consider the integral 


OO 
— OO 


OO 


z 

k=0 


f e-“' ^2n 


du 


(A. 36) 


The integral yields zero for odd values of k, thus the summation index 
can be changed and (A. 35) used to obtain 


(2m) : 

m=0 


r 2m -u^ J 

u e ^ du 



(A. 37) 


and 
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00 

J exp[-u^2/2] du = ^ 
—00 


exp[-a / 2 X 2 ] 


(A. 38) 


Differentiating n times with respect to a, and multiplying by (-1)^ 
gives the integral relation of (A. 30). 



APPENDIX B 


ALTERNATIVE METHOD FOR FINDING THE 
CONDITIONAL BIT ERROR PROBABILITY 


An alternative method exists for finding the conditional proba- 
bility of bit error for the correlation receiver using the GIN (Unique 
Waveform) model with the additional restriction of a single noise burst 

per bit time. The approach is to find the p.d.f. description for the 

k 

decision statistic, D , given that one burst has occurred. From (3.31), 


D**^|(NB=l) = I V F(T) - . k = 0,1 


(B.l) 


k I 

The p.d.f. for D | (NB=1) , i.e. , Pj^(a), is obtained by finding the cor- 
responding characteristic function, M^(u) and then taking the integral 
transform, viz. , 


oo 

— oo 


(B.2) 


where 


Mp(u) 


E{e^ 




(B.3) 


Assuming that the random variables I and T are statistically inde- 
pendent, (B.3) becomes 
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Mj^(u) = Pj(i) P^(t) dxdi . (B.4) 

i T 

Substituting (B.4) into (B.2) and interchanging the order of integra- 
tion gives 


oo 

Pd(“) ^ J' J' J"® dudxdi 


1 T 


■//'. 

i X 


(i) p^(x) 6(a - iVF(x) + C^) dxdi . (B.5) 


Using the receiver decision rule for message symbol m^ yields 


00 

P^[e^|mQ] = P[D® > 0 |nb=1] = J' Pjj(a) da 


OO 

= f fpi^i) P^(t) J' ^(a 

lx 0 


iVF(x) + Cq) dadxdi. (B.6) 


The integral of the delta function will equal one if 


iVF(x) - Cq > 0 , 


(B.7) 


otherwise the integral equals zero. It follows that 


■// 


p^(i) p^(x) dxdi 


(B.8) 


where is the region such that iVF(x) > An analogous develop- 

ment for message symbol m^ results in 
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Pl[ejnl"ii] = PtD < 0|NB=1] 




(i) p^(t) dTdi , 


(B.9) 


where is the region such that iVF(t) < C^. Since (B.8) and (B,9) 
are equivalent to (3.35), the remainder of the error analysis is un- 
changed. 



APPENDIX C 


COMPUTER PROGRAM: GINIMP 


This appendix contains the flow chart, Fortran-IV program 
listing and description for an IBM System/360 computer program 
which implements the GIN model performance analysis derived in 
Chapter 3. The program name, GINIMP, is an acronym for Generalized- 
Impulse-Noise Implementation of Model Performance. 


C.l GENERAL INFORMATION 


The performance analysis for is implemented for a speci- 

fied data bit, signal set, ENP value, noise waveform, and weighting- 
factor probability density function (p.d.f.). The user can specify 
the data bit, ENP value, and one of five weighting factor density 
functions through input data cards. The signal set and noise wave- 
form are specified by means of input data and a function subprogram, 
GRANF(x). Five typical signal sets (Unipolar, Polar, ASK, PSK, FSK) 
have been incorporated into the program. Others can be added at the 
user’s discretion by modifying the subprogram. The input data format 
and GRANF construction are explained in Section C.2. 

The normalized probability of error, from (3.42), is given by 


oo 




I . /a 
min 


(C.l) 
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Numerical techniques are used to approximate the integral solution, 
Specifically, values for the p.d.f, Pjg(n) and the time-range function 
T ^ (n) are calculated for values of n between a minimum value CXN, 
i.e., ^ maximum value XNM. An estimate to the integral is 

obtained by using the trapazoidal rule with upper and lower bound 
calculations. An explanation of how these steps are performed is 
given in Section C.2. 

C.2 FLOW CHART, PROGRAM LISTING AND EXPLANATION 

The flow chart of the GINIMP computer program is shown in Fig. 
C-1. A program listing is presented in Fig. C-2 which is found at 
the end of this Appendix. The numbers positioned beside the flow- 
chart symbols indicate the level of statement numbers used in the 
main program. 

In the first step, two data cards are read and interpreted ac- 
cording to the format and data identification given in Tables C.l 
and C.2. The input ENP and SIGSET values are used to calculate the 
FENP values , where 


FENP 


C 


0 


(l-p)Va 


(C.2) 


It has been shown in Section 3.3 that for selected signal sets this 
is a function of the ENP value. The signal-set parameter SIGSET 
determines the function to be used, cf., (3.55) and (3.56) for 
specific examples. In addition, parameters used in other parts of 


the program are calculated and initialized. 
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Fig. C-1. GINIMP Flow Chart 
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TABLE C.l 


DATA CARD 1 FORMAT AND IDENTIFICATION 


SYMBOL 

COLUMN 

LOCATION 

FORMAT 

EXPLANATION 

KBIT 

1-2 

12 

KBIT^TB = Starting time of data bit being 
analyzed 

TB 

3-14 

F12.6 

Bit duration (in ms) 

TF 

15-26 

F12.6 

Noise Waveform duration (in ms) 

ENP 

27-38 

F12.6 

Energy-to-Noise Parameter = E^VCTb/c^^ 

M 

39-40 

12 

Specifies p.d.f., Gaussian M=l, Rayleigh 
M=2, Lognormal M=3, Uniform M=4, Hyperbolic 
M=5 

PA 

41-52 

F12.6 

Parameter value associated with p.d.f. 

D 

53-64 

F12.6 

Parameter associated with noise waveform 

FREQ 

65-76 

F12.6 

Parameter associated with noise waveform 

NENP 

77-80 

14 

Number of ENP values used 


TABLE C.2 


DATA CARD 2 FORMAT AND IDENTIFICATION 


SYMBOL 

COLUMN 

LOCATION 

FORMAT 

EXPLANATION 

DF 

1-5 

F5.2 

Specifies even step increments for Slice 
Level (DF < 0) , or weighting factor 
(DF > 0) 

KKK 

6-9 

14 

Number of interpolated Fjd) values 

KK 

10-13 

14 

Number of calculated Fx(T) values 

II 

14-17 

14 

Number of Weighting-Factor values (and 
Slice Levels) used 

ISLM 

18-21 

14 

Specifies calculation of CXN (ISLM < 0) 

SLM 

22-33 

F12.6 

Maximum Slice-Level value to be used 

XNM 

34-45 

F12.6 

Maximum Weighting-Factor value to be used 

SIGSET 

46-48 

F3.1 

Signal Set to be analyzed. Unipolar = 1, 
Polar = 2, ASK = 3, PSK = 4, FSK = 5 
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The calculation of the Fj^(T) data is accomplished in the second 
step using (3.48), i.e., 

(KBIT+1)TB 

f<t-T)[sj(t)-s„(01 dt . (C.3) 

(KBIT)TB 


The program calculates KK+1 values for the FI array corresponding to 
F (T ) for the values of T =(J-1)TB/KK, where J=l,2,..., (KK+1). 

J. J J 

Since f(t-T)=0 for t<T, it follows for the first bit that KBIT=0 and 
the integral in (C.3) becomes 

TB TB 

Fi(T) = J GRANF(t)dt = 

0 


/ 


GRANF(t)dt 


(C.4) 


where 


f(t-T)[s (t)-s (t)] 
GRANF(t) = ( l-p)v ~ ~ • 


A similar result holds for (t-T)>TF. The proper limits of integration 
(A to B) are determined by subroutine LIMINT which is called using the 
statement CALL LIMINT (KBIT, TB,T,TF, A, B, ) . GRANF is provided as a func- 
tion subprogram, and for the example shown in Fig. C-2 contains the 
exponentially-decaying-sinusoid noise waveform, i.e.. 


f(t-T) = C sin[2Trf^(t-T) ] exp[-d(t-T)] (C.6) 

Input data is used to specify f^, d, and the time duration by the use 
of the mnemonics FREQ, D and TF as defined in Table C.l. The subprogram 
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calculates C such that f(t) has unit energy. The SIGSET parameter, 
defined in Table C.2, is used to determine the signal set difference 
[s^(t)-SQ (t) ] / [ (1-p) V] . The user must modify the f(t) and signal set 
descriptions if additional forms are to be used in GRANF(t). The inte- 
gration of GRANF for each T is accomplished by the IBM subroutine 

J 

QATR [44]. The accuracy of each integration is specified by FPS, and 
the subroutine returns an error message (lER) to indicate if the 
accuracy was obtained. A printout of the header I and lER is produced. 
All values of lER that indicate the accuracy FPS was not obtained are 
also printed. If all integrations are within the required accuracy 
the statement "ALL INTEGRATIONS WITHIN SPECIFIED ACCURACY" is printed. 
The maximum value stored in the FI array, FMAX, is determined and at 
the user’s option (ISLM_<0) used to calculate the maximum slice level 
SLM. 

The third step of the program, i.e., the 300 block, calculates 
the beginning n-value CXN, the step size DXN for equal steps of n 
values, and CC for equal slice-level step values. Printout includes 
headers, input data and calculated parameters. 

Values for n and the slice levels are calculated and stored in 
the XN (or SL) array in the fourth step. The user specifies by the 
value of DF whether the program uses equal step-size values for XN 
(DF<0) or equal step-size values for SL (DF^O) . 

During the fifth step the time-range function T - (n) is approxi- 

R1 

mated by calculating values for the TR array. This is accomplished 
by comparing each value in the FI array to each value in the SL array 
and increasing the number in the TR array by one if F1(I)>SL(J). A 



95 


linear interpolation is performed between the calculated FI values 
and each point is classified as above. The elements of the TR array 
can be divided by the number of FI points considered to yield an 
approximation to Tj^^(n)/T^. 

The probability of error is calculated in the last operation. 
First, the TR array is modified to obtain T^^(n)/T^, and the result 
multiplied by the corresponding p.d.f. value obtained from the func- 
tion subprogram PDF(X,M,PA). The user selects one of five p.d.f. 
descriptions by selecting the proper m-parameter value as listed in 
Table C.l. The PA value is a parameter used by either the lognormal 
(PA>0) or the hyperbolic 3<PA<6 distributions. The numerical- 
integration subroutine QTFGB is then used to calculate the approxi- 
mation of Pt [ e ]. The subroutine is similar to the IBM subroutine 
1 m 

QTFG [44], but has been modified to include upper and lower bounds 
for the integral. 

A printout is made from the various arrays, the most signifi- 
cant of which are the cumulative probability of errors obtained 
from the integration subroutine, namely POEL the lower bound, POE 
the probability of error, and POEU the upper bound. The last value 
in the POE array represents checks for accuracy can 

be performed. First, the values of POE should converge to the last 
value. If this does not happen, the user should increase the maximum 
value of n, XNM, and repeat the computation until convergence is 
obtained. Second, the lower and upper bounds can be compared to the 
POE value to determine if a sufficient number of XN values (II) were 


used. 
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DIMENSION TR(500) ,SL(500), ANP(500) tAUX(500) , FI (500) 
DOUBLE PRECISION CXN , DXN , XN ( 500 ) , POE ( 500 ) , POE L ( 500 ) 
DOUBLE PRECISION P0EU(500) 

EXTERNAL GRANF 

COMMON T,TF,FREQ»D, ISIG 

READ(5,101)KBIT,TB,TF,ENP ,M , PA ,D »FREO,NENP 

101 FORMAT! 12 , 3F 1 2 . 6, 1 2 ,3F 12 .6 , 1 A ) 

READI5,102)DF ,KKK,KK, 1 1 1 I SLMf SLM ,XNM , SI GSET 

102 F0RMAT(F5.2,4I4,2F12.6,F3.1) 

IRE=0 

111=500 
ISIG=SIGSET 
DO 605 IENP=1,NENP 
IF ( lENP.EQ. 1 ) GO TO 107 
ENP=10.*ENP 
107 CONTINUE 

GO T0(103tl0A,105»103, 103), ISIG 

103 FENP=SQRT(ENP/2. ) 

GO TO 106 

104 FENP=SQRT(ENP) 

GO TO 106 

105 FENP=SQRT(ENP/4. ) 

106 CONTINUE 

IF ( IENP.GE.2 ) GO TO 301 

C SETS UP CONDITIONS, DETERMINES LIMITS, CALCULATES F(T) 

C FOR EACH VALUE OF T SPECIFIED, AND STORES RESULTS IN 

C THE FI ARRAY 

C 

FMAX=0. 

DTT=TB/(1000.«KK*KKK) 

KK 1=KK+1 
KKK1=KKK+1 
WRITE (6,200 ) 

200 FORMAT ( '0',4X,' I',4X,'IER',/, • «#***=!:*****#«* • , 

1 • S|c :;=:{! :(S5;C5;: #5j::{c *5):# :{c #5): sic# !!:#:!£ 5jt!{: ^:{:j!c<c:!c5!:s!: , ) 

DO 204 1=1, KKl 
T= ( I-l )*TB/KK 

CALL LIMINT(KBIT,TB,T,TF,A,B ) 

FPS=1.0E-7 
A=A*.001 
B = B=!=.001 

CALL QATR( A,B,F PS , I I I , GRAN F , Y , I ER , AUX ) 

FI ( I)=Y 

IF(IER) 201,203,201 

201 CONTINUE 
IRE = IRE + 1 

WRITE (6,202 ) I , lER 


Fig. C-2. GINIMP Program Listing 



202 F0RMAT(2X,2I4) 

203 CONTINUE 
YABS=ABS (Y) 

IF( YABS.LT.FMAX )G0 TO 204 
FMAX=YABS 

204 CONTINUE 

IF(IRE) 508»506t508 

506 WRITE(6,507) 

507 F0RMAT(2Xt'ALL INTEGRATIONS WITHIN SPECIFIED ACCURACY') 

508 CONTINUE 

IF ( ISLM) 205 ,205,301 

205 SLM=FMAX+. 000001 

C CALCULATES PARAMETERS BASED ON INPUT DATA, 

C HEADER AND DATA OUTPUT 

C 

301 CXN=FENP/SLM 
DXN=(XNM-CXN )/( II-l) 

TR ( 1)=0. 

SL( 1 )=0. 

XN ( 1 1 ) = 1 0=«‘*4 
XN( 1 )=CXN 
SL ( II )=FENP/CXN 
CC=SLM/(II-1) 

WRITE(6,302) SIGSET 

302 FORMAT! '!• ,20X, 'PERFORMANCE ANALYSIS FOR GIN MODEL *, 

1 ' AND SIGSET = ' ,F3.1,/,2X, ' , 

1 !###«#», /,2X, 'THE INPUT PARAMETERS ARE',/) 

WRITE(6,304) 

304 F0RMAT(6X, ' KBIT' ,8X, 'TB' , lOX , ' TF ' , 9X , ' ENP ' ,9X,'M',9X, 
l'PA',10X, 'D' , lOX, 'FREQ' ) 

WRITE(6,305 ) KBIT , TB , TF , E NP , M, PA ,D , FREQ 

305 FORMAT (7X, 12 ,5 X , 3F 12 . 8 , 2X, I 2 , IX , 3F 12. 4 , / ) 

WRITE (6,306 ) 

306 FORMAT (9X, 'DF ' , 8X, ' KKK ' , IX , ' KK ' ,2X , ' I I ',4X, 'CXN' ,9X, 

1 'DXN' ,9X, 'SLM' ) 

WRITE (6,307)DF, KKK,KK, I I, CXN, DXN, SLM 

307 F0RMAT(6X,F 12.8,3I4,3F12.8) 

WRITE(6,309) 

309 F0RMAT(2X,12(/ ),2X, 'NOISE WAVEFORM USED*«*«' , 

I'EXPONENTIALLY DECAYING SINUSOID') 

C CALCULATES EVEN SL OR EVEN XN 

C 

IF(DF )403 ,401,401 
401 DO 402 J=2, II 
TR( J)=0. 


Fig. C-2 (Continued). GINIMP Program Listing 
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402 

403 


404 

405 

c 

c 


501 

502 


503 

504 

505 
C 4 

c 

c 

601 


604 


603 

605 


SL( J ) = ( J-1 )*CC 

XN ( II+1-J) = FENP/SL( J ) 

GO TO 405 
DO 404 I=2»II 
TR( I )=0. 

XN ( I)=XN (I-l )+DXN 
SL( I I + l-I )=FENP/XN( I ) 

CONTINUE 

INTERPOLATES AND CLASSIFIES THE F(T) VALUES 


DO 505 I=2»KK1 

IF( I.LT.KKl) GO TO 501 

JJK=KKK1 

GO TO 502 

JJK=KKK 

DO 504 J=1,JJK 
TT=(J-1)#DTT 

F2=ABS( (F1(I)-F1(I-1))#( (J-1 )/KKK)+Fl (I~l) ) 

DO 504 L = l, II 

IF(F2.GE.SL(L) ) GO TO 503 

GO TO 504 

TR(L )=TR( L )+l 

CONTINUE 

CONTINUE 

sjsjicjic :Jc sis 3^:^ jJ; sissCssicsicsississis^sIssisstcsjcsIssicsic 

CALCULATES THE PROBABILITY OF ERROR 


WRITE(6,601) 

FORMAT! • 1S80X»'CUMULATIVE PROBABILITY OF ERROR*, /»5X, 
1 ' I* , lOX, • XN* , 14X, 'SL* , 14X, *TR' ,13X, *ANP' , 12X, • LOWER* , 


1 IIX, * AVERAGE* , 9X , ' UPP ER * , / , 2X , • *=i=***«**4:^:*:!:***«#^=*## » 

«JU >V ft 

V ^1% ^|V >1^ >1^ V 


t 

t 


■ o 


ft \ 

> ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ • j 


DO 604 1=1,11 

TR ( I ) =TR { I ) / (KK*KKK + 1 ) 

INEW=II+1-I 

YPDF = PDF (XN ( INEW ) ,M,PA) 

ANP( INEW)=TR( I ) *YPDF 
CONTINUE 

CALL QTFGB(XN,ANP,POE,POEL,POEU,II ) 

WRITE(6,603 ) (I ,XN( I ),SL( I I + l-I ),TR( II + l-I) ,ANP( I ), 
IPOEL ( I ) ,POE( I ) , POEU( I ) , 1 = 1 ,I I ) 

F0RMAT(2X,I4 ,7F16.9) 

CONTINUE 

STOP 

END 


Fig. C-2 (Continued). GINIMP Program Listing 
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SUBROUTINE QT FG B { X , Y , Z » ZL , ZU ,ND I M ) 

DIMENSION Y 1500) 

DOUBLE PRECISION X ( 500 ) , Z ( 500 ) » Z L ( 500 ) , ZU ( 5 00) 

DOUBLE PRECISION SUM2 , SUM L2 , SUMU2 

SUM2=0. 

SUML2=0. 

SUMU2=0. 

IF (NDIM-1)4,3,1 
1 DO 2 I=2»NDIM 

SUM1=SUM2 
SUML1=SUML2 
SUMU1=SUMU2 
DELTAX = X{ I )“X ( I-l ) 

SUMP=.5=i=DELTAX*(Y( I )+Y( I-l) ) 

SUM2=SUM2+SUMP 
DELTAY=Y ( I )-Y( I-l ) 

DELY2=Y( I+l)-Y( I ) 

DELYO=Y( I-l )-Y(I-2 ) 

DELYY2 = DELY2-DE LTAY 
DELYY1=DELTAY-DELY0 
S1=DELY0/{X(I-1 )-X(I-2) ) 

S2=DELY2/(X ( I+l)-X( I ) ) 

IF(S1.EQ.S2) GO TO 10 

IF { I.EQ.2.0R.I .EQ.MDIM) GO TO 8 

D = (DELTAY-S2*DE LTAX ) / ( S 1-S 2 ) 

DI =DELTAX-D 

A=( Y( I-l)+D*Sl/2. )*D+( Y( I)-DI*S2/2. )*DI 
IF (DELYY2.GE.0..AND.DELYY1.GE.0. ) GO TO 5 
IF(DELYY2.LT.O. .AND. DELYY1 .lt. 0. ) GO TO 6 
SUML2=SUML2 + SUMP 
SUMU2=SUMU2+SUMP 
GO TO 7 

5 SUML2=SUML2+A 
SUMU2=SUMU2+SUMP 
GO TO 7 

6 SUML2=SUML2+SUMP 
SUMU2=SUMU2+A 

GO TO 7 

8 Al=( 1-2 )=^ ( Y(NDI M-1 )+Sl*DELTAX/2. )/(NDIM-2) + 

1 (I-NDIM)*(Y{2)-S2*DELTAX/2. )/(2-NDIM) 

A=A1*DELTAX 

IF (DELYY2.GE .O..AND. I .EQ. 2) GO TO 9 

IF(DELYY1 .GE.O. .AND. I.EQ.NDIM) GO TO 9 

SUML2=SUML2+SUMP 

SUMU2=SUMU2+A 

GO TO 7 

9 SUML2=SUML2+A 
SUMU2 = SUMU2+SUMP 


Fig. C-2 (Continued). GINIMP Program Listing 
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GO TO 7 

10 SUML2=SUML2+SUMP 
SUMU2=SUMU2+SUMP 
7 ZL( I-l )=SUML1 
ZU ( I“1) = SUMU1 

2 Zn-1)=SUM1 

3 Z(NDIM)=SUM2 
ZL(NDIM)=SUML2 
ZU (ND IM ) =SUMU2 

4 RETURN 
END 

« 5^ j!« !(: sjc j!s # Jjc # 4c * Jjt sis # sic ### sjc sis 5): sj: Jjc ## sj: ## J}: # 5}!* «#*#*#* ****** * 

SUBROUTINE L IM I NT ( KB IT , TB» T , TF , A, B ) 

J=TF/TB 

IFIKBIT. GT.O. .AND. KBIT. LT.J) GO TO 50 
IF (KBIT.EQ.O..AND.J.EQ.O. ) GO TO 53 
IF(KBIT.GT.O.. AND. J.EQ. KBIT) GO TO 55 
IF (KBIT.EO. ( J+1. ) ) GO TO 51 
IF(KBIT.EQ.O..AND.J.GT.O.) GO TO 57 
GO TO 58 

50 A=KBIT#TB 

B= (KBIT+1. )*TB 
GO TO 58 

51 IF ( (T+TF ).LE.KBIT*TB ) GO TO 52 

56 A=KBIT*TB 
B=T+TF 

GO TO 58 

52 A=0. 

B=0. 

GO TO 5 8 

53 IF( (T+TF) .LE.TB ) GO TO 54 

57 A = T 
B=TB 

GO TO 58 

54 A=T 
B=T+TF 
GO TO 58 

55 IF {(T+TF ).LE.( (KBIT+1.)=!=TB) ) GO TO 56 
GO TO 50 

58 RETURN 
END 


Fig. C-2 (Continued). GINIMP Program Listing 
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FUNCTION GRANF(X) 

COMMON T,TF,FREQ,D,ISIG 

F=FREQ 

TFF=TF*.001 

PI =3.14159265358979 

A=l./( PI#F )«SORT( (D=!'(D«#2+4.#(PI*F)*«2 ))/( 1-EXP (- 
12.*D4TFF))) 

GRANF = A*SIN(2.*PI=f'F^(X-T’!'.001) ) «EXP {-D*( X-T*.001 ) ) 

GO T0(l,l,2t2t3),ISIG 

1 RETURN 

2 GRANF=GRANF#SIN(2000.*PI«X) 

RETURN 

3 GRANF=GRANF«(SIN(2000.=!'PI*X)-SIN(4000.=«'PI=s=X) ) 

RETURN 

END 

Q*#*#:{ts}:3}cj}:s!5:{£ sic 45 4es}c4!5!e:}ts}c4c si: ❖ 

C 

FUNCTION PDF(X,M,PA) 

PI=3. 14159265358979 
GO T0( It2»3t4,5 ) tM 

1 PDF = ( 1/SQRT (2.*PI ) )#EXP(-X«=S'2/2. ) 

RETURN 

2 PDF=ABS(X)=»=EXP(-X»'«2) 

RETURN 

3 PD1 = 1/(SQRT (8.#PI )#PA#ABS(X) ) 

PD2=EXP(-( ALOG( ABS ( X ) ) +PA**2 ) «#2/ ( 2.*P A**2 ) ) 

PDF=PD1=«=PD2 

RETURN 

4 Z=ABS(X) 

IF(Z-SQRT(3.) )41,41,42 

41 PDF = 1/SQRT( 12. ) 

RETURN 

42 PDF=0 
RETURN 

5 PD3 = ( PA-1 )=«= ( (PA-2. )=«=(PA-3. ) )*#( (PA-1 )/2. ) 

PD4=(SQRT( (PA-2 . )*(PA-3. ) )+SQRT(2. )*ABS(X) )=s=*PA 
PDF=1/SQRT( 2.)*PD3/PD4 

RETURN 

END 


Fig. C— 2 (Continued). GINIMP Program Listing 



APPENDIX D 


COMPUTER PROGRAM: LIMITER 


The Fortran-IV computer program LIMITER implements the perfor- 
mance results of Section 5.2. The program name is a mnemonic to 
emphasize the association with the generalized limiter of Fig. 5.2. 
A succinct explanation of this program is presented by the use of a 
flow chart and program listing. 


D. 1 GENERAL INFORMATION 


The program implements (5.12) to calculate P^[e^|m^] for a 
given data bit, signal set, ENP value, noise waveform, limiter 
parameters, and weighting factor p.d.f. The use of either positive 
or negative values for the weighting factor must also be controlled. 
These parameters are specified by input data and function subprograms. 
Four typical signal sets (Unipolar, Polar, ASK and PSK) and the five 
weighting-factor p.d.f.s of Table 3.1 have been incorporated into the 
program. Numerical techniques are used to approximate the integral 
solutions. The minimum weighting-factor value can be specified or 
calculated in the program by an iteration process. An explanation 
of the program operation is given in Section D.2. 
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D.2 FLOW CHART, PROGRAM LISTING AND EXPLANATION 

The LIMITER program flow chart is shown in Fig, D-1, A program 
listing is presented in Fig, D-2, which is located at the end of this 
Appendix. The first symbol represents the operations performed by 
statement numbers 0 through 299. Five data cards are read by the 
format and data identification given in Tables D.l through D,5, Two 
arrays are calculated based on the function subprogram specification 
of the signal set and noise waveform. The ENP value is selected from 
the ENPA array and the program is initialized for positive or negative 
weighting-factor calculations. The calculation of Cq/(cV) is made 
according to the signal set selected. Headers and data are printed 
and parameters are initialized for the weighting factor iteration 
process . 


TABLE D.l 


DATA CARD 1 FORMAT AND IDENTIFICATION 


SYMBOL 

COLUMN 

LOCATION 

FORMAT 

EXPLANATION 

ENPA(l) 
ENPA(2) 
ENPA (3) 

2-13 

14-25 

26-37 

F12 . 10 
F12 . 10 
F12 , 10 

ENPA is a one-dimensional array with 
three elements containing the ENP 
Values 


The second step in the flow chart corresponds to operations per- 
formed by statements with numbers 300 through 899. Normally, the 
iteration process for the minimiim positive (or maximum negative) 
weighting factor will be utilized (by specifying ISCH = 1) . The 
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time-range function is calculated for a weighting-factor value of 
CXO/2. Successive approximation steps in the weighting-factor value 
are made in the iteration attempt to locate the largest value that 
gives a zero value for the time range. After this operation is com- 
plete, the array of weighting-factor values is finalized. The non- 
linear operation {Nf(t-T)}°^^ is performed in the series of steps 
through the 500-level statement numbers. These results are mul- 
tiplied by the signal-set array data and integrated in the 600-level 
statements. An approximation to the time-range function is made in 
the 700-level statements by classifying and interpolating the integral 
results that exceed Cq/CoV). 


TABLE D.2 


DATA CARD 2 TORMAT AND IDENTIFICATION 


SYMBOL 

COLUMN 

LOCATION 

FORMAT 

EXPLANATION 

SIGSET 

3-13 

F10.4 

Signal Set to be analyzed. Unipolar =1, 
Polar = 2, ASK = 3, PSK = 4 

FI 

14-23 

F10.4 

Frequency of s;L(t) for ASK and PSK 

FO 

24-33 

F10.4 

Not used, intended for frequency of 
SgCt) for FSK. 

KBIT 

34-37 

14 

KBIT*TB = Starting time of data bit 
being analyzed 

MPDF 

38-41 

14 

Specifies p.d.f., Gaussian = 1, 
Rayleigh = 2, Lognormal = 3, Uniform = 
4, and Hyperbolic = 5 

PA 

42-51 

F10.8 

Parameter value associated with p.d.f. 

TB 

52-61 

F10.8 

Bit duration (in sec) 


Multiplication of the time-range data by the corresponding p.d.f. 


values is accomplished by the 900-level program statements. Headers 
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and output data are printed. An estimate of the error probability 
is calculated in the last step. The results are presented graphically 
with the final value in the graph representing the estimate. Checks 
for convergence of the integral and correct weighting-factor step size 
should be made. 


TABLE D.3 


DATA CARD 3 FORMAT AND IDENTIFICATION 


SYMBOL 

COLUMN 

LOCATION 

FORMAT 

EXPLANATION 

NIFF 

3-6 

14 

Number of interpolated points between 
each time-shift point 

NTL 

7-10 

14 

Number of sample points used in bit time 

Tb 

NTLT 

11-14 

14 

Number of NTL samples per time shift 

NXN 

15-18 

14 

Number of Weighting-Factor values used 
between minimum and maximum values 

CONN 

19-28 

F10.3 

Defines LIMITER operation, 
3n = /CONN Un 

CONP 

29-38 

F10.3 

Defines LIMITER operation. 
Bp = /CONP Up 

XNMP 

39-48 

F10.4 

Maximum positive Weighting-Factor value 

XNMN 

49-58 

F10.4 

Minimum negative Weighting-Factor value 
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TABLE D.4 


DATA CARD 4 FORMAT AND IDENTIFICATION 


SYMBOL 

COLUMN 

LOCATION 

FORMAT 

EXPLANATION 

CXNP(l) 

CXNN(l) 

CXNP(2) 

CXNN(2) 

CXNP(3) 

CXNN(3) 


IQQQH 

CXNP is a one-dimensional area with three 
elements containing the minimum positive 
Weighting-Factor values for each of the 
three ENP values. CXNN is a similar array 
for the maximum negative Weighting-Factor 
values . 


TABLE D.5 


DATA CARD 5 FORMAT AND IDENTIFICATION 


SYMBOL 

COLUMN 

LOCATION 

FORMAT 

EXPLANATION 

MJ 


14 

Message symbol to be analysed mQ=0,m2=l 

ICXN 


14 

Specifies sign of Weighting-Factor values 
Positive = 1, Negative = 2 

ISCH 

11-14 

14 

Specifies use CXNP or CXNN (ISCH=0) or 
iteration routine to calculate value 
(ISCH=1) 

CXO 

15-26 

F12.6 

Maximum value of iteration step 
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6 

1 

2 

5 

4 


3 


C 


10 1 


20 2 


203 

204 


DI i'-'ENSIOi'l XN ( 600 ) ,T IME (6 0 0) »XNF ( 1200 ) ,SXMF (6 00) , 

ISA (600 ) ,YY (101 ) ,Z (600) , ZZ( 100) , Ai'lS ( 100 ) , ENP A ( 3 ) » CXMiM ( 3 ) , 
IFA ( 1200) ,CXi'!P(3) 

COlH-iOfJ F1,F0 

PI =3, 14159265358979 

P^EAD (5 ,6) ( ENPA ( I ) , 1 = 1, 3 ) 

FORMAK 1X,3F12.10) 

READ (5 , 1 ) S IGSET , F 1 , FO , K B IT , MPDF , P A , TR 
FORMAK 2X,3F10.4, 214, 2F10 .8 ) 

R E A D ( 5 , 2 ) N I P F , N T L , N T L T , N X M , C Q i\i N , C G N P , X M W P , X M M M 
FO RM A T ( 2 X , 4 I 4 , 4F 1 0 . 3 ) 

READ ( 5 , 5 ) ( CXNP ( I ) , C XNM ( I ) , I = 1 , 3 ) 

FORMAK 6F10 .7) 

READ (5 ,4)MJ, ICXN, ISCH,CXD 
FORMAK 2 X, 3 I 4, FI 2. 6 ) 

uDTTC/A.ci\CTr:CE:T.c:i .i 


0RMAK2X,3 I4,F12.6 ) 

RITE (6,3 ) SIGSET,F1,F0,KBIT,MPDF,PA,TB, 
XM ,CONN ,CONP , XNMP , XNMM 

r“ r'l D f .'i A "T / f 1 I o V I *■'- so V*- 

I U K r I A I \ * J. ’ J ^ ^ ^ *,•» -o 

^ ^ O# 0 > ^ 1 ^ 


f^lIPF ,MTL,MTLT, 


# / I /A \/ 

^ 1 -.. ^ 1 %. -r|- ^ ^ ^ 1 N A y 


1 ' PERFORMANCE ANALYSIS USING L I M I TER ' , / , 3X , ' S I G SET ' , 


14X, »FK , 


8X, *F0',8X, 'KBIT ' ,1X, 'MPDF ' , IX 


, 'PA ' ,9X, 


'TB 


1/,2X,3F10.4,2I4 ,2X,F10.8,1X,F10.8, //,2X, 'NIPF' , IX, 
1 'NTL' ,1X , 'NTLT ' , IX, 'NXN' , 4X , 'CONN' ,6X, 'CONP' , 6X , 
I'XNHP' ,6X, 'XriMN ' , / , 2X , 4 14, 4F 1 0 . 3 ) 

ISIG=SIGSET 
DTL = TB/(NTL-1 ) 

CALCUATE SA AND FA ARRAYS 

DO 101 1 = 1, NTL 

TL = ( I-l ) =:=DTL 

TIME ( I )=TL 

SA ( I)=S(TL, ISIG) 

TTL= (KBIT-D-TB+TL 

TTTL=TTL+TB 

FA( I )=F (TTL) 

FA (NTL+I )=F ( TTTL ) 

SELECT ENP, CXN, XNM AND FENP 
PRINT HEADER AND INPUT DATA 
INITIALIZE EOR CXN ITERATION 
DO 910 I ENP =1,3 
IF( ICXN-1 )202,202,203 
CXN=CXNP ( lENP ) 

Xi\'M = X(liiP 
GO TO 204 
CXN=CXNN( lENH ) 

Xi\'M=XNMN 
CONT INUE 
ENP=ENPA ( lENP ) 

GO TO (210,211, 212,213,214 ), ISIG 


Fig. D-2. LIMITER Program Listing 
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210 FENP = S0RT ( ENP/2 . ) 

GO TO 215 

211 FENP=SQRT ( ENP ) 

GO TO 215 

212 FENP = SQRT ( ENP/4 . ) 

GO TO 215 

213 FENP = SORT (ENP/2 . ) 

GO TO 215 

214 FENP = SORT ( ENP/2 . ) 

215 CONTINUE 

WRITE(6,292) EM P , M J , ICXN ,C XN 

292 FORMAT! ' 1 'THE FOLLOWING ANALYSIS IS FOR ENP=',F12.3, 
15X, 'MJ = SI2,5X, • ICXN=' , I2,5X, ' CXN = • , F 1 2 .6 , / ) 

CXS=CXO 
I CYXN=0 
ISTEP=0 
IXFX=ISCH 

C ITERATION FOR CXN (PART 1) 

868 CXS=CXS/2. 

ISTEP=ISTEP+1 

WRITE ( 6,872) I STEP, I CYXNt CXN 

872 F0RMAT(2X,'THE',I3,2X,'STEP YIELDS I CYX N = • , I 5 , 2X , ' AND ' , 
1' CXN=',F12.6) 

IF ( ISTEP.EQ.50) GO TO 910 

IF( ISTEP.EQ.2.AND.ICYXN.E0.0) GO TO 910 

IF ( ICYXN >866,866,867 

866 CXN=CXN+CXS 
GO TO 860 

867 CXN=CXN-CXS 

860 DXN=( XNM-CXN ) / (NXN-9 ) 

C CALCUATE XN VALUES 

DO 201 IXN=1,NXM 
I CYXN=0 

IF ( IXN.GE.l .AND. IXN.LT.ll ) GO TO 250 
XXN=( IXN-9)*DXN+CXN 
GO TO 251 

250 XXN=( IXN-1 )=!=(DXN/9. )+CXN 

251 XN(IXN)=XXN 

C MULTIPLY EACH FA BY XN 

DO 301 1=1, NTL 
XNF ( I ) =XXN*FA( I ) 

INTL=NTL+I 

301 XNF ( INTL)=XXN*F A( INTO 

C TIME SHIFT XN*FA 

NT=(NTL/NTLT)+1 
DO 401 J=1,NT 
IT=( J-D-NTLT 

C THE INDEX IT CORRESPONDS TO THE TIME SHIFT T 
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DO 405 L=1,NTL 
405 SXNF (L )=XNF(NTL-IT+L) 

C CLIP SHIFTED XN#FA 

DO 509 L = 1»NTL 
TL = (KBIT-1)*TB + (L-1) =^-DTL 
IF(MJ-1 )507t508 »508 

507 GO TO {5071 ,5072,5073,5074,5075) tlSIG 

5071 DUML=0. 

DUMU=SQRT (2 .-ENP/ ( TB**2 ) ) 

GO TO 510 

5072 DUML=0. 

DUMU=SQRT (4.*ENP/(TB=<'*2 ) ) 

GO TO 510 

5073 DUML =-SQRT (4.*ENP/ (TB*«2 ) ) 

DUMU=SQRT(4 .*ENP/ (TB«*2 ) ) 

GO TO 510 

50 74 DUML = -SQRT( 2.4ENP/ (TB«*2 ) ) * (1 .-S I N ( 2 .*P I *F 1*T L ) ) 
DUMU=SQRT ( 2.*ENP/ (TB**2 ) ) * (1 .+S I N ( 2. * P I *F 1* TL ) ) 

GO TO 510 
5075 DUML=0. 

DUMU=0. 

GO TO 510 

508 GO TO (5081 ,5082,5083,5084,5085) ,ISIG 

5081 DUML=-SQRT(2.*ENP/(TB**2) ) 

DUMU=0. 

GO TO 510 

5082 DUML=-SQRT(4.*ENP/(TB**2) ) 

DUMU=0. 

GO TO 510 

5083 DUML=-SQRT (4.*ENP/ ( TB**2 ) )* (1 .+S I N ( 2 . *P I *F 1 *TL ) ) 
DUMU = S0RT(4 .*ENIP/ (TB**2) ) * (1 . -S I N ( 2 . *P I «F 1 *TL ) ) 

GO TO 510 

5084 DUML=-SQRT(2.«ENP/(TB**2) ) * (1 . +S I N ( 2 .*P I *F 1«T L ) ) 
DUMU=SORT( 2.-ENP/(TB*=:'2 ) )’- ( 1 .-S IN ( 2.*P I*F 1«TL ) ) 

GO TO 510 

5085 DUML=0. 

DUMU=0. 

510 CONTINUE 

IF (SXNF ( D.GE .DUML.AND.SXNFd ) .LE.DUMU) GO TO 517 
IF(SXNF (L ) .GT.DUMU.AND.MJ.EO.O) GO TO 513 
IF (SXNF(L).GT.DUMU.AND.MJ.EO.l) GO TO 514 
IF(SXNF (L) .LT.DUML.AND.MJ.EO.O) GO TO 515 
IF (SXNF(L).LT.DUML.AND.MJ.EQ.l) GO TO 516 
513 GOTO (5131,5132,5133,5134,5135 ), ISIG 

5131 SXNF(L)=S0RT(2.*C0NP4ENP/ (TB**2) ) 

GO TO 517 

5132 SXNF(L )=SQRT (ENP/ (TB«*2 ) ) *( SORT(CONP )+l . ) 

GO TO 517 
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5133 SXNF (L )=S0RT(4.*C0NP*ENP/(TB*#2) ) 

GO TO 517 

5134 SXNF (L )=SQRT(2.*ENP/(TB*#2 ) ) * ( SORT ( CONP ) + 
1SIN(2.-PI*F l^^TL) ) 

GO TO 517 

5135 SXNF(L)=0. 

GO TO 517 

514 GO TO (5141 ,5142t5143, 5144, 5145) tlSIG 

5141 SXNF (L)=SQRT(2.*ENP/(TB**2 ) ) * ( SORT { CONP )-l . ) 
GO TO 517 

5142 SXNF (L)=SQRT(ENP/(TB*42 ) )* ( SORT ( CONP )- 1 . ) 

GO TO 517 

5143 SXNF (L )=S0RT(4.«ENP/(TB*«2 ) ) * ( SORT ( CONP ) - 
1 SIN (2.*P I«F 1*TL) ) 

GO TO 517 

5144 SXNF(L) = S0RT(2.=«'ENP/ (TB*#2) )*(SORT(CONP )- 
1SIN(2.*PI*F14TL ) ) 

GO TO 517 

5145 SXNF(L)=0. 

GO TO 517 

515 GO TO (5151,5152,5153,5154,5155 ), ISIG 

5151 SXNF(L)=-SQRT(2.*CONN«ENP/(TB#*2) ) 

GO TO 517 

5152 SXNF(L) = SQRT(ENP/(TB**2) ) ^'( -SORT ( CONN ) + 1 . ) 

GO TO 517 

5153 SXNF(L)=-SQRT(4.=!'C0NN*ENP/(TB#*2) ) 

GO TO 517 

5154 SXNF(L)=SORT(2.*ENP/(TB**2) )*( -SORT ( CONN ) + 

IS IN ( 2.-PI^'Fl*TL ) ) 

GO TO 517 

5155 SXNF(L)=0. 

GO TO 517 

516 GO TO (5161,5162,5163,5164,5165), ISIG 

5161 SXNF(L)=S0RT(2.-ENP/ (TB**2) )*( -SORT ( CONN )- 1 . ) 
GO TO 517 

5162 SXNF(L)=S0RT(ENP/(TB**2) )*( -SORT ( CONN )- 1 . ) 

GO TO 517 

5163 SXNF(L)=S0RT(4.*ENP/(TB**2) )*(-SORT(CONN)- 
1SIN(2.*PI*F1*TL ) ) 

GO TO 517 

5164 SXNF (L )=S0RT(2.-ENP/(TB<'*2 ) )*( -SORT ( CONN ) - 
1SIN(2.*P I*F 1*TL) ) 

GO TO 517 

5165 SXNF(L)=0. 

517 CONTINUE 
509 CONTINUE 

C SXNF NOW REPRESENTS SHIFTED AND CLIPPED XNF 

C 
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C MULTIPLY CLIPPED XN*FA BY SA 

C INTEGRATE TO OBTAIN FJ(N,T) AND STORE IN YY ARRAY 

DO 601 M=1,NTL 

601 SXNF (M)=SXNF(M) «SA(M) 

CALL QTFG(T IME ,SXNF»Z,NTL ) 

IF(MJ)602»602,603 

602 YY(J)=Z(NTL) 

GO TO 401 

603 YY ( J) =-Z (NTL ) 

401 CONTINUE 

C INTERPOLATE AND CLASSIFY FJ(N,T) VALUES TO OBTAIN 

C TIME RANGE DATA 

DO 701 IC=2tNT 

IF(YY( IC-1 ).LT.FENP.AND.YY( lO.LT.FENP ) GO TO 701 
IF ( IC-NT )704 ,705 ,705 

705 NIP=NIPF 

C NIPF IS THE NUMBER OF INTERPOLATED POINTS FIXED 

GO TO 706 
704 NIP=NIPF-1 

706 DO 703 ICC=1,NI P 

C NIP IS NUMBER INTERPOLATED POINTS 

YDUM=( (YY( IC)-YY{ IC-1) )*{ICC-1) ) / ( N I P F -1 ) +Y Y ( I C- 1 ) 

IF (YDUM-FENP )703,707,707 

707 I CYXN= ICYXN+1 
703 CONTINUE 

701 CONTINUE 

C ITERATION FOR CXN (PART 2) 

YK=ICYXN 

IF ( IXFX.EQ.O)GO TO 970 

IF( ICYXN.GT.O.AND.ICYXN.LE .4)G0 TO 871 
GO TO 868 
871 IXFX=0 

C • MULTIPLY TIME RANGE DATA BY PDF VALUES 
C PRINT HEADER AND OUTPUT DATA 

970 ANS(IXN)=(YK/( (NIPF-1)=MNT-1)+1) )«PDF(XN( IXN) ,MPDF,PA) 

IF( IXN-1)973,973,974 

973 WRITE(6,975) 

975 FORMAT (3X, ' IXN‘ , 15 X , ' YK ' , 1 8X , • ANS ' ) 

974 WRITE{6,991) I XN , YK , ANS ( I XN ) 

991 FORMAT (2X, 14 , 10 X , F 1 2 .6 , lOX , F 1 2. 6 ) 

201 CONTINUE 

C CALCUATE PROBABILITY OF ERROR 

900 CALL OTFG(XN,ANS,ZZ,NXN) 

CALL GRAPHK XN, ZZ,NXN, 2 ) 

WRITE(6,499) ( YY ( J ) , J =1 , N T ) 

499 FORMAT (2X, lOFll .6) 

910 CONTINUE 

STOP 
END 
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FUNCTION S(TL,ISIG) 

COMMON Fl,FO 

PI=3. 14159265358979 

GO TO ( 1 , 2,3 ,4,5) , ISIG 

1 S=l. 

RETURN 

2 S=1 . 

RETURN 

3 S=SIN(2.=:^PI*F1*TL) 

RETURN 

4 S=SIN ( 2.-=P I-F1*TL ) 

RETURN 

5 S=SIN( 2.-PI*Fl=:<TL)-SIN (2.*PI*FO*TL) 
RETURN 

END 




^ ^ 


FUNCTION F(X) 

IF(X)9,8,8 

8 CONTINUE 
TFF=.004 

IF (X-TFF )11,11,9 
11 D=500. 

C FREO=G 

G=1000. 

PI =3.14159265358979 

A = l,/{ PI*G I'-^SORT ( (D*(D**2+4.=;'(PI«G)**2 ))/( 1 -EXP ( -2*D=!'TFF ) ) ) 
F = A='^SIN ( 2.*PI*G*X )*EXP (-D*X ) 

GO TO 10 

9 F = 0. 

10 RETURN 
END 
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